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Abstract 

The Lagrange reconstructing polynomial [Shu C. W.: SIAM Rev. 51 (2009) 82-126] of a function fix) on a given set 
of equidistant (Ajc = const) points {x, + £Ax; I e {-M_, ■ ■ ■ , +M+}} is defined as the polynomial whose sliding (with 
x) averages on [x - ^Ax,x + \Ax] are equal to the Lagrange interpolating polynomial of f(x) on the same stencil 
[Gerolymos G.A.: J. Approx. Theory 163 (2011) 267-305]. We first study the fundamental functions of Lagrange 
reconstruction, show that these polynomials have only real and distinct roots, which are never located at the cell- 
interfaces (half-points) x, + rijAx (n e Z), and obtain several identities. Using these identities, we show that there 
exists a unique representation of the Lagrange reconstructing polynomial on {/ - M_, ■ ■ ■ ,i + M+} as a combination of 
the Lagrange reconstructing polynomials on Neville substencils [Carlini E. , Ferretti R. , Russo G. : SIAM J. Sci. Comp. 
11 (2005) 1071-1091], with weights which are rational functions of ^ (x = x, -H^Ax) [Liu Y.Y., Shu C.W., Zhang M.P.: 
Acta Math. Appl. Sinica 25 (2009) 503-538], and give an analytical recursive expression of the weight-functions. 
We show that all of the poles of the rational weight-functions are real, and that there can be no poles at half-points. 
We then use the analytical expression of the weight-functions, combined with the factorization of the fundamental 
functions of Lagrange reconstruction, to obtain a formal proof of convexity (positivity of the weight-functions) in the 
neighborhood of ^ = 5' iff of the substencils contain either point / or point / + 1 (or both). 

Keywords: reconstruction, (Lagrangian) interpolation and reconstruction, hyperbolic PDEs, finite differences, finite 
volumes 
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1. Introduction 

Polynomial interpolation and/or polynomial reconstruction are the basic numerical approximation operations in- 
volved in the development of WENO schemes iUISIj which are widely used ^ for the discretization of (hyperbolic) 
PDEs, particularly when the solution contains discontinuities. Following Godunov's theorem [4], these schemes intro- 
duce nonlinearity in the approximation (with respect to the reconstructed function h{x) or to its cell-averages fix)), 
to combine high-order with monotonicity. Central to the development of these methods ijsl [ll is the underlying lin- 
ear approximation, where the interpolating and/or the reconstructing polynomial on a given stencil is 
represented by a combination of the corresponding (interpolating or reconstructing) polynomials on substencils. We 
introduce the following definitions 

Definition 1.1 (Stencil 10, Definition 4. 1, p. 283]). Consider a 1-D homogeneous computational mesh 

Xi = xi -H (/ - l)Ax Ax = const e M>o (la) 

Assume 

M:=M_+M+>0 (lb) 
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The set of contiguous points 

S/,M_,M+ := {i - M_, ■ ■ ■ , / + M+} (Ic) 

is defined as the discretization-stencil in the neighborhood of /, with M_ neighbors to the left and M+ neighbors to the 
right. The stencil S, m_,m+ (flcl l contains M + 1 > points and has a length of M intervals. If M-t > then the stencil 
contains the pivot-point /. If M_M+ < then the stencil does not contain the pivot-point We will note 

[Si,M-,Mj ■^[Xi-M-,Xi+Mj C M (Id) 

the interval defined by the edge-points of the stencil. □ 

Definition 1.2 (Neville substencils). Let S,-,m_,m+ be a discretization stencil on a homogeneous grid (Definition II. Il l 
with 

M M_ H- M+ > 2 (2a) 

Assume 

K,<M-l ; K.eNo (2b) 

The -H 1 > 1 substencils 

Si,M^-k,,M,-K,+h -.^{i-M^+k,,--- ,i + M+-K,+k,} y k,e{0,--- ,K,} (2c) 
each of which contains M - + 1 points and which satisfy 

Si,M_-k„M^-K,+k, - ^i,M-,M+ (2d) 

{.^nis <^ < • ,^ V4,Ws e {0, ■ ■ ■ .iiTsI (2e) 

S,-M_+fc+i,/+M+-A:,+fc+i = (s,-M_+fc,/+M+-A:,+fc \ {i - M_ H- ^j) U {/ + M+ - K, + k, + 1} Vk, e {0, ■ ■ ■ , K, - I] (2f) 
are the (M-K, + l)-orde£| substencils of S,;m_,m+, corresponding to the /Tj-level subdivision of S,-,m-,m+- D 



Definition 1.3 (Reconstruction pair |l7|, Definition 2.1, p. 270]). Assume that Ax e M>o is a constant length, and that 



the functions / : / — > M and h : I — > M are defined on the interval / = [a - 4 Ax, H- 4 Ax] c K, satisfying everywhere 



/(x) = — /2(^)fl'^ Vx6[fl,/7] (3a) 

Ax J.i-i Ax 

assuming the existence of the integral in ( [3a] i. We will note the functions /(x) and /!(x) related by ( l3at 



=-R(I;A.)(/) (3b) 
/ (3c) 

and will call / and h a reconstruction pair on [a, b], in view of the computation of the l-derivative|l □ 



'in the sense that the Lagrange interpolating and reconstructing polynomials on each of the substencils )2ct are 0(Ax'^ *^**')-accurate approx- 
imations [3, Proposition 4.6, p. 289]. 

+ -Ax) - h^"^'\x - -Ax) 

2 By (1 Lemma 2.2, p. 271], (3al =^ f-"\x) = V.v e [a,/)] Vn e |1, ■ ■ ■ ,^1, exactly, assuming f(x) and 

Ax 

h(x) are of class C'^[a — jAx, b + jAx\. 
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Definition 1.4 (Lagrange reconstructing polynomial iQ Definition 2.3, p. 271]). Let pi,M_,M+(xi + f Ax; x,, Ax; /) de- 
note the Lagrange interpolating polynomial [8, pp. 186-189] of the real function / : M — > M on the stencil S,,m_,m+ 
(Definition ! 1. 11 1. Its reconstruction pair (Definition ll.3t 



PR,.M-.Mjxi + ^Ax; X/, Ax; /) [^(1;a.v)(P/,a/_,m+)](-«; + ^Ax; x,-, Ax; /) (4) 

will be called the Lagrange reconstructing polynomial on the stencil S, m_,m+ ■ D 

We study representations where the polynomial approximation on S,;m_,m+ (Definition 11.11 ) is expressed as a 
weighted sum of the corresponding polynomial approximations on the /Tj + 1 substencils (Definition ! 1.21 ) 

PR^,M-MMi + ^^'^^i^^x'J) ^^CTR^^M^^M^^K^kSS) PRi,M_-k,M^-KMMi + ^Ax; X;, Ax; /) (5a) 

Pi,M-,MMi + ^^^'^ Ax; /) = ^ o-i^M-^M^^K„kM) PiM--k,M^-K,^kMi + ^Ax; X,-, Ax; /) (5b) 

with weight-functions {o-iM_.M^.K^,kSS) 'he interpolation case or o-R^^M_M+,K„kSS) the reconstruction case) which 
are independent of the approximated function (/(x) in the interpolation case or h{x) in the reconstruction case). The 
subscripts M±, and ks in (|5]l indicate that the weight-functions depend on the stencil (M±), on the level of subdi- 
vision {Ks) and on the particular substencil (^s)- Because the weight-functions are independent of the approximated 
function (/(x) or h{x)) they are usually called linear weights |3]. Alternatively, since the weights combine the inter- 
polating (or reconstructing) polynomials on the substencils to exactly the interpolating (or reconstructing) polynomial 
on the entire stencil, they recover the highest possible accuracy (between weighted combinations of the substencils) 
and, for this reason, they are alternatively called optimal (in the sense of accuracy) weights 121 |9l]. 

The underlying linear interpolation or reconstruction used in WENO |5, 3] schemes on the general stencil {/ - 
M_,--- ,/ + M+) (Definition !!. 1! ) can be obtained by writing the approximation error fli (56a), p. 292] for the 
/Ts + 1 (^s 6 {0, ■ ■ ■ , ^Ts)) substencils {/ - M_ + k^, - ■ ■ ,i + M+ - + h] (Definition !!.!! ). each of which has an 
error of (9(Ax''^"*-'^') |7, Proposition 4.7, p. 292]. At any fixed point x, -H ^Ax, we can in this way construct a 



{K^ + 1) X (Ks + 1) linear system {eg HlOl (13), p. 8489]) for the weights which linearly combine the approximated 
values on the substencils to obtain an (9(Ax*'^')-accurate approximation, recovering the accuracy (and indeed the 
exact value |5, 3]) of the entire stencil {/ - M_, ■ • • , / -i- M+), at the chosen fixed point x, + ^Ax. It is known by 
numerical experiment Isllil, that, for stencils symmetric around x,- (/e M_ = M+) these linear or optimal weights 
for the {K^ - M_ = M+)-subdivision, can be calculated at the fixed point ^ - \ {eg |10, Tab. 3, p. 8484]), ie for 
this choice of {M+,Ks,^] the linear system [10., (13), p. 8489] is not singular Shu [5] has given examples of other 
choices of {M+, K^,^) for which the linear system is singular Obviously the weights are functions of f , parametrized 
by [M^,KM 

The Neville-Aitken algorithm |8, pp. 204-209] constructs the interpolating polynomial on {/ - M_, ■ ■ ■ ,i + M+}, 
by recursive combination of the interpolating polynomials on substencils, with weights which are also polynomials 
of X 13, PP- 204-209]. Carlini et al. |6], working on the Lagrange interpolating polynomial in the context of centered 
(central) WENO schemes, recognized the connexion between the Neville-Aitken algorithm |^ pp. 204-209] and 
the determination of the optimal weights, and gave the explicit expression [6, (3.6,4.10), pp. 1074-1079] of the 
polynomial weight-functions o"/ r_i,,y _i_i,(^) which combine the Lagrange interpolating polynomials on the Ks + I - 
{r—\)+\ substencils {/- (r- 1) -H^s, •■•,' + '"-('"- 1) + ^s} to obtain the Lagrange interpolating polynomial on the big 
stencil {; - (r - 1), ■ • ■ ,i + r} which contains an odd number of M = 2r - 1 intervals and an even number of M -H 1 -2r 
points. This result was also confirmed by Liu et al. (2.2), p. 506] who further gave the analytical expression ifTll 
(2. 1 8), p. 511] for the polynomial weight-functions cri r,r,r,kM) which combine the Lagrange interpolating polynomials 
on the Ks + l = r+l substencils {i-r + k^, - ■ ■ J + r-r + kf^} to obtain the Lagrange interpolating polynomial on the big 
stencil {/-r, ■ ■ • ,i + r] which contains an even number of M -2r intervals and an odd number of M+1 = 2r+ 1 points. 
For both cases it is shown [4. JJJ that e [-1,1] the linear weights are positive, and as a consequence the above 



combination of substencils is convex € [-1, 1]. In a recent work !|12|| we extended these results for the general 



A's-level subdivision of an arbitrary stencil X,_m-,/+m+ '■- [xi-m-, ■ • • ^-'Ci+M+l c R of M + 1 := M_ -i- M+ + 1 distinct 
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ordered points on an inhomogeneous grid to A's + 1 < M substencils ^{-m^+kj+m^-k.+k '■- {^i-M_+k,, ■ ■ ■ , Xi+M. -K. +L] 
(ks e {0, ■ • ■ , Ks]), and used a general recurrence relation [12, (4e), Lemma 2. 1] to obtain a simplified expression lll2l 
Proposition 3. 1] of the weight-functions for the Lagrange interpolating polynomial ( fSbl i. and to prove positivity in the 
interval x e [jc,-m-+a's-i , -«/+m+-a:s+i] which contains at least 1 cell (at least 2 grid-points) iff < lfl2L Proposition 
3.2]. 

Looking more carefully into (l5bl l one notices that it is directly related to Miihlbach's theorem |fl3l Theorem 2. 1 , p. 
100], corresponding to lll3i (2.2,2.3), p. 100]. Miihlbach [13] expresses the coefficient cri M_,M+,K^,kSO (l5b] i in terms of 
quotients of determinants of interpolation-error functions, directly obtained by the Cramer solution [14, Proposition 
5.1.1, p. 72] of error-eliminating linear systems. Miihlbach 1 13] studies Chebyshev-systems satisfying interpolatory 
conditions. In the reconstructing polynomial case (ISal l. the usual linear system approach (13), p. 8489] is 
equivalent to the algorithm of Miihlbach [13, Theorem 2.1, p. 100] with the important difference that in (l5al l we study 
polynomials Pr^(x) whose linear functionals Pijx) - [R^-^.^^^{pi{^)](x) (Definition ! 1.3b satisfy interpolatory conditions, 
so that the existence and uniqueness proofs in |13l Theorem 2.1, p. 100] are not directly applicable. Nonetheless, the 
general recurrence relation for weight-functions proven in llT2l (4e), Lemma 2.1], only requires that the (Ks - l)-level 
subdivision can be defined. Therefore, finding a general expression for the weight-functions ctr^ m-.m^.k^mS^) < l5a] l. is 
tantamount to solving the problem of the {K^ - l)-level subdivision for the Lagrange reconstructing polynomial. 

Although the reconstructing polynomial Isl IH [tI] is even more widely used in WENO discretizations, the develop- 
ment of practical WENO schemes [QHIi^l^ invariably followed the aforementioned linear system approach [10, (13), 
p. 8489], using symbolic calculation. There is little analytical work on the weight-functions cr^,,M_,M+,A:„<rs(f) which 
combine the Lagrange reconstructing polynomials on the K^ + l < M substencils {/ - M_ +ks,--- ,i + M+ -K^+ k^] to 
obtain the Lagrange reconstructing polynomial on the big stencil (/ - M_, • • ■ ,i + M+V Only recently, Liu et al. Ill | 
studied particular families of stencils and subdivisions, using symbolic computation^ Liu et al. |ll[l have concen- 
trated on the usual WENO substencils0 In the reconstruction case, it was shown by construction [Tljthat the optimal 
weight-functions are not polynomials, as in the interpolation case [|6, n, 12 1, but, instead, rational functions of ^ 
(x = Xj + ^Ax), implying that at the poles of these rational functions the weight-functions cannot be defined. For 
upwind-biased schemes the big stencil (/ - (r - 1), ■ • ■ , / -H (r - 1)} (r € N>2) which is centered around the 

point /, and upwind-biased with respect to the cell-face / + j, is subdivided [11, Tab. 3.2, p. 516] into + I - r 
substencils (/ - (r - X) + k^, - ■ ■ ,i + k^, A:s e {0, ■ • ■ , K^. For centered schemes, the big stencil {/ - (r - 1), ■ ■ ■ , / -H r) 
which is centered with respect to the cell-face / + | ||5lll5l|3||, and as a consequence downwind-biased with respect to 
the point /, is subdivided into /Tj -H 1 = r + 1 substencils {i -(r -X^ + k^, - • ■ ,i + k^\, A:s e (0, ■ ■ • , K^\. In [1 k (3.2-3.4), 



p. 514] an algorithm is sketched for computing the rational weight-functions, which are tabulated up to r = 7 01 1 



Tab. 3.2, p. 516] for the upwind-biased case (even number of intervals) and up to r = 6 [1 1, Tab. 3.5, p. 518] for the 
centered case (odd number of intervals). We remarked in [7, p. 298] that both these families can be grouped together 
as the subdivision of the general stencil {/ - [y J, • ■ • ,i + M- [yj} into A's H- 1 = |^yj + 1 substencils, in the range 
M e {2 ■ ■ ■ , 11}. These weights were further analyzed to determine the regions of convexity of the representation 
(positivity of the weight-functions). These important results [1 LI include explicit expressions of the weight-functions 
for the particular stencils which were studied, but a general analytical expression of the optimal weight-functions 
for the representation of the Lagrange reconstructing polynomial by combination of substencils is not yet available. 



contrary to the inte rpola ting polynomial case [|6l IHJ, [12|| . The work of Liu et al. [11] is based on the reconstruction via 



primitive approach 0161 pp. 243-244], as developed in , where the integral (primitive) h{^) d( of the function 
h(x), which is reconstructed from its sliding averages f(x) (Definition ! 1.3b . is used. 

Despite the enormous successes of the reconstruction via primitive approach fldi, pp. 243-244] in designing 
and analyzing practical WENO schemes llsl S [TTtl the reconstruction via deconvolution approach ifm 244-246] is 
conceptually more straightforward, since it directly uses the unknown function which is reconstructed from cell- 
averages, and sometimes simplifies analytical work. In a recent work [7, Lemma 2.5, p. 272] we have provided the 
analytical solution of the deconvolution problem !fl6[ (3.13b), p. 244], which expresses the unknown function h(x). 



^ Liu et al. [11] have examined, using symbolic calculation, the computation and positivity of linear (optimal) weight-functions in WENO 
interpolation, reconstruction and integration. 

"^In the nomenclature of Shu used in Liu et al. ITlll . stencils are defined in tenns of cell-interfaces (half-points), and the term nodes in 

fllll denotes cells, so that the stencil {i - (r - 1), ■ ■ ■ ,i + (r- 1)| is defined in 1 11, Tab. 3.2, p. 516] as [i — r + j, ■ ■ ■ ,i + r — and the stencil 
{/- (r- 1), ••■,! + fl is defined in |TT, Tab. 3.5, p. 518] as {i - r + i, ■ ■ ■ , i + r + ^ | 
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which is reconstructed from its sliding averages, as a series of the derivatives of the sliding averages ^ (10b), 

p. 272]. This analytical solution of the deconvolution problem 1 16, (3. 13b), p. 244] allows the analytical computation 
of the approximation error of the Lagrange reconstructing polynomial |7, Proposition 4.7, p. 292], which would have 
been necessary to build the general linear system ifiol (13), p. 8489] for the weight-functions o^R^^M-.M^.K^,kS^) < l5a] i. 

In the present worlj| we use relations and concepts developed in |7], along with the general recurrence relation 
for the generation of weight-functions proven in L12i, Lemma 2.1], to extend the analysis of Liu et al. |11], both 
by providing general analytical expressions (and existence and uniqueness proofs) of the rational weight-functions, 
but also by studying the general case of the subdivision of an arbitrarily biased stencil on a homogeneous grid, 
{i - M_, ■ ■ ■ , / + M+) (M± E Z : M :- + > 2) containing M intervals, into + I < M substencils 
(Definition ll.2l ) of equal length of M - intervals, each shifted by 1 cell with respect to its neighbors {K^ is free to 
take all possible values £ {1, ■ ■ ■ ,M- 1)). We also prove several relations concerning the Lagrange reconstructing 
polynomial. 

In ij2]we very briefly summarize those results for the Lagrange reconstructing polynomial and its approximation 
error obtained in ['t'] which are necessary in the present workO 
In m we study the basis polynomials oii^ M-,M+ 

/(f) {{ £ {-M_,--- ,M+) which L7» Proposition 4.5, p. 287] 
represent the Lagrange reconstructing polynomial on S,;m_,m+ := {' - M^, ■ ■ ■ ,i + M+), with coordinates the values 
fi+i :- f(xi + {!Ax) of the cell-averages of the reconstructed function. These results, which include an analysis 
of the roots of aR^^M_^M^/(^) and relations with the polynomials /lR,,M_,M+,«(f) appearing in the expression of the 
approximation error of the Lagrange reconstructing polynomial 10, Proposition 4.7, p. 292] are the starting point 
for the construction of the Lagrange reconstructing polynomial as a combination of the Lagrange reconstructing 
polynomials on substencils. 



In fj4]we use the results of ^to establish a 1 -level subdivision rule (Lemma [4.2l i. by which, applying 11121 Lemma 
2.1], we construct (Proposition l4.5t an analytical recursive expression of the weight-functions for a general subdivision 
of an arbitrarily biased stencil on a homogeneous grid. We prove the uniqueness of the rational weight-functions 



(Proposition I4.7l i. and we show by studying their poles (all of which are real) that it is always possible to define 
the weight-functions at half-nodes - n + j , n e Z). Finally, we prove (Theorem 14. 14) the convexity of the 
representation (l5a] i in the neighborhood of f = 5, for all subdivisions (Definition 1 1.2l i for which all of the substencils 
contain either point ; or point / + 1 (or both). 

2. Reconstruction background 

In a recent work |7] we have studied the exact and approximate reconstruction of a function h(x). We have 
obtained the general analytical solution of the deconvolution of Taylor-series problem [16. (3.13), pp. 244-254], and 
used this solution in developing analytical relations for the approximation error of polynomial reconstruction on an 
arbitrary stencil in a homogeneous grid fl\. We briefly summarize those results of [T] which are the starting point of 
the analysis presented in the present work, and which are necessary for completeness. 

Lemma 2.1 (Derivatives of reconstruction pairs). Let h - R(\;hx){f) be a reconstruction pair (Definition 11.31 ), and 
assume that f(x) and h{x) are of class C^[a - ^Ajc, b + |Ax]. Then 

h = R(v,Uf) ^ h^"^ = /?(i;Ax)(r) Vn £ {1, • ■ • , A^) (6) 
Proof. We have by direct integration 

h{x + IrAx) - h(x - \Ax) J /m. 
Ax Jr-lAr Ax 



1 px+ ^ Ax 



from the fundamental property of reconstruction pairs [0 Lemma 2.2, p. 271], proving (|6]l for n = 1, and by induction 
Vn E {I,--- ,A^). □ 



5 In §6.1, pp. 297-300], we had sketched, without giving any proof or analysis, some of the problems which are solved in the present paper 
Furthennore, at that time, we had not proven the conjectured convexity i3.[3l lTlll . in the neighborhood of f = j. 

In the present work we also make extensive use of relations concerning reconstruction pairs Definition 2.1, p. 270], and in particular 
polynomial reconstruction pairs Theorem 5.1, p. 296]. 



Expressions for the Lagrange interpolating polynomial and its approximation error are widely available in the 
literature fs", pp. 186-189], and they are only included in the following to define notation, for completeness, but also to 
highlight analogies and differences between Lagrange interpolation and Lagrange reconstruction. In [7, Propositions 
4.5, 4.6, 4.7] we developed corresponding analytical expressions for the Lagrange reconstructing polynomial, which 
can be summarized as 

Proposition 2.2 (Lagrange polynomial interpolation and reconstruction on S,jw _m+ ll3l)- M± e 7, : M :- 

M_ + M+ > ( llbl ). Let h — R(i-Ax)(f) be a reconstruction pair (Definition 11.31 ). Then the Lagrange interpolat- 
ing (pi,m_mJ ^f^d reconstructing (pri,m_,m^) polynomials of f(x) (Definition 1 1 .41 ) on S,_m_,m+ (Definition II . II ) are 

(.X; +^Ax;x;, Ax;/) = ^ q'r,,m_,m^/(^) /(jc; + Mx) (8a) 

A't, 

^h(xi+^Ax)+ 2 ;U«,,M_,M,,„(a Ax«/(«>(x,) + 0(Ax^"+i) (8b) 

n=M+l 

^h(xi+^Ax)+ 2 iR„M_,M„„(a Ax"/!<">(x, +^Ax) + 0(Ax^-+i) (8c) 

n=M+l 



Pi,M-MMi + ^Ax; Xi, Ax; /) = ^ cei,M-M^AO f(xi + (Ax) (9a) 

e=-M- 

A't, 

^f(xi+^Ax)+ 2 /i/,M_,M,,„(a Ax" /"^U) + 0(Ax^"+') (9b) 

A't, 

=/(x, +^Ax)+ 2 4m_,a/„„(^)Ax«/("'(x, +^Ax) + 0(A/'"+1) (9c) 



«=M+1 



and, provided that h(x) is sufficiently smooth Vx e [x,_m_ - jAx, x,+jvf+ + |Ax], their approximation errors are de- 



fined by dSBISbll, or equivalently by (l8c]|9c]). The functions OR^^M^^M^^ei^), fiR,,M-,M^,n(^)' ^.M-.m+A^I 0!i,m_,m^/(^), 
l-ii,M-,M+,n(0 '^^d ^i,M-,M+,ii(0 are polynomials in 



X - Xi 



Ax 



with coefficients depending only on M± and are defined by 



M 



m=() \ k=0 

LfJ 



T2k(m + 2ky M -I 

; Im )m+2k+\.e+M- + l 

ml 



M 



(10) 



-M_</'<+M+ (lla) 



fiRi,M-,M- 



k=Q 
n-M-\ 



m=a\ k=0 



s\ ml 



s > M + 1 



(lib) 

(11c) 
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ai,M-M+A^) --^(m!^ )m+l/+M_+l 



m=0 

1 

S 



m=0 



M^<e < +M+ 



s>M+\ 



f=o ^• 



n>M+\ 



(12a) 
(12b) 
(12c) 



where (^*V ^ are the elements of the inverse Vandermonde matrix on S,;m_,m+ IlL Definition 4.3, p. 283], expressed 



by 10, (43a,43b), pp. 283-284]]] VM_,M+,m,s are defined 



^M-,M+,m,k • 



m+I,f+M- + I 



A:€ No 



(13) 



a«c/ f/je numbers t„ satisfy^ 

T0=1 



T2k 



-T2.S 



, , 22-^ (2i + 



^ 22*-2-^ {2k -2s+ 1)! 22-^ (2s + 1)! 



T2n+1 = 



k>0 
n>0 



(14a) 

(14b) 
□ 



Remark 2.3 (Alternative expressions for Q;/^M_,M+,f(^) ( I12al ) and a'R,,M_,M+/(^) ( II lab ). The polynomials ai,M_,M+AO 
the fundamental functions [Sl pp. 183-197] of Lagrange interpolation on the stencil S,;m_,m+, and can also be ex- 
pressed, upon replacing Xf - x,- + £Ax in (9.4), p. 184], as 



k=-M_ 



(15a) 



' By H Lemma 4.4, pp. 283-284] 



/t- 1 
(• - 1 



By R Lemma 4.4, pp. 283-284] 



V(,;e{l,--- ,M+1) 



0<k<M 
0<m< M 



Vfe€No 



' By 0, Theorem 2.9, pp. 275-276] the numbers r„ can be defined as t„ := — ^?r"'(0) from the generating function grix) '■ 



sinh 5 j: 
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Shu isl (2.19), p. 336] has shown, using the reconstruction via primitive approach lfl6l pp. 243-244], that the 
fundamental functions ( II 1 al l of the Lag range reconstructing polynomial on S,,m_,m+, PRi,m_,m+(-^i + ^Ax;jic,, Ax;/) 
:an equivalently be expressed ao 

M+l M+l 

p=0 q=0 

BM+1 P*!" q*m 
(2.19)] V 



=f+M_ + l 



]~[(m-/5) 



p=0 



□ 



Remark 2.4 (Mapping /?(i;Ax))- By |0, Theorem 5.1, p. 296] the mapping ^(1;Aa) is a bijection of the (M + 1)- 
dimensional space Mm[-'c] of polynomials of degree < M onto itself. This implies that polynomial reconstruction 
pairs are unique li Lemma 3.1, p. 296]. Let p{x) e Rm[x] be a polynomial of degree < M. Then by |7, Theorem 5.1, 
p. 296] 

( q(x) [/?(i;A.)(/?)](x) € Rm[x] 
V/7(x)eRMW ^ ] deg(^) = deg(p) (16) 

[ coefr[x'''=sW,;?(jc)] = coefr[jc'''=gW, ^(x)] 

where the property that the linear operator 7?(i;A.i) '■ Rm{x\ — > Rm[jc] conserves the degree of the polynomial follows 
directly from [T* Lemma 3.1, p. 277]. Furthermore, the coefficients of the leading power of p{x) e Mm[-«L and of 
q(x) :- \R(i:A.x-)(p)](x) £ Km[.';] can be shown to be equal, by straightforward application of the expression JtI (26f), 
Lemma 3.1, p. 277] One consequence of these properties is that several relations obtained for the interpolating 
polynomial have their direct analogues for the reconstructing polynomial and vice-versa. □ 



3. Fundamental polynomials of Lagrange interpolation and reconstruction 

The construction of a recursive formulation ( |j4|i for the linear weight-functions (l5al i is based on the repre sentations 
of the Lagrange reconstructing polynomial (l8at and of its approximation error (l8ci . It is therefore necessary to gain 
some insight on the fundamental functions ( II lal ) of Lagrange reconstruction, and on the truncation-error polynomials 
(fTTcl i. 

3.1. Reconstruction pairs of fundamental polynomials 

Each of the fundamental polynomials of Lagrange reconstruction aR^ M-.M^AO ( II lab is intimately related to the 
corresponding fundamental polynomial of Lagrange interpolation ai^M_,M+AO (I12al ). as can be seen by using ( |9al[8al l 
in the reconstruction-pair-defining relation (l3at . 

Proposition 3.1 (Reconstruction pairs aR^^M-,M+,t - ^(1;1)(Q'/,m_,m+/))- Assume M± e Z : M := M_ -i- M+ > ( fTbl l. 
The polynomial a/t^^M-.M+.d^^ dl lab appearing in the representation dSab of the reconstructing polynomial on the sten- 
cil S,-,M_,M+ (Definition ll.lb is the reconstruction pai^^of the corresponding polynomial ai^M^,M+,eiO I ll2ab appearing 



The coiTespondence of the present indicia] notation with the one used by Shu f?, (2.19), p. 336] is rshu = ishu = ^ + M_, <:shii = M + 1 = 
M- + M+ + 1, and we use again x = Xi + ^Ax ilQ\ . 

" By Lemma 3.1, p. 277] if deg(p) = M, and p{x, + ( Ax) = Eff-o'V,,, then q(x, + ( Ax) := [R(i;Ax)(p)](-r, + ( Ax) = Zff-oCg™ i'" and 
by 0, (26f), p. 277] 

[1(260] 1 ^ /„^^M.D^ 

"l'^ = M! ^ ^2*' W2t (^ + 2*:)! - c„„ 

A:=0 

^^on a unit-spacing grid, Ax = I 



in the representation ( |9al ) of the interpolating polynomial on the same stencil 

J (--2 



V {e {-M_,--- ,M+ 



(17) 



Proof. By Definition 1 1.41 of the reconstructing polynomial, we have, using (Bal 

-X, +f A.r+ J A.V 



1 r- 

Ax J.v, + 



PRxM-mA^'^ Xi, Ax; /) t/^ (18a) 

and using the representation (l8at of the reconstructing polynomial and the representation (l9at of the interpolating 
polynomial, we readily obtain by (llSal i 



/(x,. + ^Ax)™GM0 



V^€ R 
Vx/ e M 
VAx € R>o 
V/:R^ 



(18b) 



□ 



Since ( llSbl l is valid V/ : R — > R, it proves, using the definitions (l3al|3b]l, (fTTT l. 

Lemma 3.2 (aR^M-.M+.K^) Om„[^] aj^M-,M+AO)- Assume M± e Z : M := M_ + M+ > ([lb]l anc/ let 

M 

m=0 

denote the 0-element of the space RM[f] of all polynomials of degree < M. None of the polynomials ( II lal I12al ), of 
degree M in ^ is identically 



,M. 



(19b) 



Furthermore 



coefr[^*^,aR,,M_,M,/(f)] = coefr[^'^',a/,M_,M,/(^)] = (-1) 



deg[a«,,M_,M+/(^)] = deg[Q'/,M_,M+/(f)] = M 



■M 



M 



Ml \r + M_ 



(19c) 

Vfe{-M_,--- ,M+) (19d) 



Proof. It is well known, and also obvious from the expression (115a) . that the fundamental polynomials oj^m-.m+A^) 
ofLagrange interpolation on a stencil of M+1 equidistant points are OM„[f](f). Since bv (ITtTi m m, A^) + Or„[^](^) 
is equal to the definite integral of or^^m-M+A^) the interval - 5, ^ + Al, it follows that or^^m-.m^A^) 
By ( II lal I12ab we have ( |19c| ). It is easy to show by direct computatioq^j that coeff[^*^, m_,m+,;(^)] is given by 



By Proposition l2.2l 



COeff[f*^,aKj,M_,M„f(f)] 



1^ M-M J 



Q ... ..M+l-M-l 



I (43a)l 



-(0^^ )jM+l,f+M_ + l 



7,(43b)l + l 



2j a-- 



/t- 1 



{k- \y.\e + M- + \ - I 



k- 1 
M+ 1 - 1 



/ M 1 


M 




M 



1 / M 



■ (-1)' 



1 / M 



#0 Vfe{-M- 



< M 



where we used the expressions [7, (43a,43b), pp. 283-284] for the elements of the inverse of the Vandermonde matrix (fr^] P-EJ- well known 



properties of the unsigned Stirling numbers of the first kind (17, Tab. 264, p. 264], m < n ^0 - 



V m, n e N and 



1 V ;i e No. 



( |19d| ). By ( fT6b this is also the coefficient of of the polynomial ai,M_,M+,e(0 ( I12al ). which (Proposition 13. lb is the 
reconstraction pair of Q;R,,M_,M+,f(^), proving ( |19(H . □ 



• ,M+)}). Assume Mj, e Z : M := M_+M+ > GB. The{M+\) 
6 {-M_, ■ ■ ■ ,M+), constitute a basis of the (M + V)- dimensional 



Proposition 3.3 (Basis {q'r,,m_,m+/(^), ^ £ (-^-, 
unique polynomials Or^^M-M+AO ^ OkmK] Eil' 
space Ma/[^] o/a// polynomials of degree < M. 

Proof. This can be proved either by Proposition 13. ll or directly. 

Proof bv Proposition l3.1l It is a well-known fact iQl, and also obvious from ( I15al l. that the (M + 1) unique 
polynomials {oi^m-M+AO^^ ^ {-M_,--- ,M+)) are linearly independent and span the (M + l)-dimensional space 
I^m[^] of all polynomials of degree < M in ^. By Lemma [TTI each polynomial aRi.M-.M+.f - R(i;Ax)(^i,M-.MtA 
( fTTb . and by OTl Theorem 5.1, p. 296] the mapping R(\-,&x) ■ ^m[^] — > is a bijection. Hence the image of 

{a/,M_,M„f(f),^ e {-M_, ■ ■ ■ ,M+)), {aR^.M-.M,A^)J ^ i'M-, ■ ■ ■ ,M+]} is also a basis of Rm[^]. 

Direct proof: Existence of the polynomials ffs, ,m_,m+/(^) satisfying dSal i was proved by construction ^ Proposition 
4.5, p. 287] yielding ( II lab . Recall that by |7, Theorem 5.1, p. 296] the linear operator R(i-ax) (Definition II. 3b is a 
bijection of the vector space Km[j^] of all polynomials of degree < M onto itself. Obviously, by ff. Lemma 3.1, p. 
277] the same properties apply to the inverse operator Rj^.^^y Since, V/7(x) € Km[ji;] p^'''\x) = > M + 1, the 
reconstructing polynomial (Definition 1 1.4b of p{x) on the stencil S,;m_,m+ {i - M^, ■ ■ ■ ,i + M+] (Definition 1 1.1b is 
exactly equal to the reconstruction pair of p{x) (Definition ll.4l) . q(x) :- [R(\-Ax){p)](x), by JHcb . By (IHal l we have 



q(x) 



Ax) 



(q)](xi + 1 Ax) 



Since (I20al) holds Vx, e 



. and VAjc e IR>o we may set x/ 



V<?(x) e Mm[^] 
Vx e M 
Vx; e M 
VAx 6 R>o 



(20a) 



and Ax = 1 in (I20ab to obtain 



= ^ aRxM-M^A^) [^(1;1)(?)](^) 



t=-M_ 



Vq(x) e Rm[x] 
VxeR 



(20b) 



By ( I20bb . the M + 1 polynomials Q'R,,Af_,M+/(-«), ^ e {-M_, ■ ■ ■ ,M+] span Rm[^], and since dim(RM[x]) = M + 1 
they form a basis of Rm[x]. They are therefore linearly independent llisll . and as a consequence OR„[f] ( I19bl ). a fact 
already proven in Lemma IT2l □ 



3.2. Roots of fundamental polynomials 

Because of ( fTTI l for every value returned by the polynomial a!M_,M+A^i) at point ^/ e R, there exists a nearby 
point 6 R such that aR^^M-,M+A^R\) - ctiM-M^ASi)^ the distance between the 2 points being - ^/| < |. This 
can be formalized as 

Lemma 3.4 (Q'/,M_,M„f([^i,6]) £ a'R,,M_,M,/((^i - + 5)) V^i,^2 e R : ^1 < ^2)- Assume e Z : M -.^ + 
M+ > (HB. r/ien 



6 {-M_, ■ ■ ■ ,M+) V^/ E R 3 ^s, € (^/ - i + i) c R : ,M_,M„f (^r, ) = a/,M_,M„f(^/) 



(21a) 



where Ur^^M-M+AO <I1 lab ant/ o:i^M-,m+AO <ll2ab ore the fundamental polynomials of Lagrange reconstruction and 
interpolation, respectively (Proposition 12 . 2b , implying that 



■it e {-M_, ■ ■ • ,M+} a/,M_,M,/([^i,^2]) £ a«,,M_,M„f((^i - ^,^2 + |)) V^i,^ e M : ^1 < ^2 (21b) 
Proof. The proof follows immediately from Proposition |TT] By (fTTb 



44 



(22) 



Using the mean value theorem for the definite integral lll9l p. 352] in (l22l) yields (12 lab , from which (I21bb is easily 
proved by contradiction. □ 
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The fundamental polynomials of the Lagrange interpolating polynomial a!.M-,M^AO < ll2al i are polynomials of 
degree M in ^ ( I12al ). and it is well known [Q] and obvious from their expression ( I15al l that their M roots are the integer 
nodes {-M_, ■■ • ,M+]\\t] 

ai,M-M.An)^0 VnE{-M_,--- ,M+)\{f} € {-M_, • ■ • , M+) (23a) 

ai,M-.M.M)*Q V^€M\{{-M_,--- ,M+)\K}} V^e{-M_,--- ,M+) (23b) 

The fundamental polynomials of the Lagrange reconstructing polynomial or^^m^m+A^) (11 lalllSbT i are also polynomi- 
als of degree M in ^ ( lllal i. but the expressions ( I Hal llSbl i are too complicated to directly give information about their 
roots. It is nonetheless easy, using Lemma [T4l to show that 



Proposition 3.5 (Roots of aR^^M-,M+,A^) dHab ). Assume M± eZ:M:=M_+M+>0 ( fTbl i. The M roots of the 
degree M in ^ polynomials Q'r,,m_,m+/(^) (II lab are all distinct and real, and there is exactly 1 root in each open 
interval (n - i n + i) V n e {-M_, ■ ■ ■ , M+} \ {{}, ie 



Vn6{-M_,--- ,M+)\K) 



3! ^R,,M-,M^,e,n e(n- ^,n + ^) ciR : aR,,M_,M+,f(f«,,M_,M^_/,n) = (24a) 



a«,,M_,M,,f(^) ^0 6 M \ {^R,,M_,M„f,„; n e {-M_, ■ • ■ , M+) \ {^)| e {-M_, • ■ ■ , M+) (24b) 



Proof. The proof follows immediately from Lemma l34l by writing (12 lab at each of the M roots ( I23ab of aiM_,M+,ciO- 
By Proposition 13 .21 the polynomial Q's,,m_,m+/(^) ?^ OR„[f](^). Furthermore deg[Q'R,_M_,M+/(^)] = A# ( I19cb . and since 

there are exactly M := M_ + M+ elements in |{-M_, ■ ■ ■ , M+} \ {^}|, the roots ( I24ab are, by the fundamental theorem 

of algebra and its corollaries lHoi pp. 282-289], the only roots of afi,,M_,M+/(^), which proves (I24bb . and uniqueness 
(3!) in ( I24ab . by contradiction. □ 



Remark 3.6 (Extrema of OR^M^MtAS) ( II lab ). It is straightforward to show that each fundamental polynomial of La- 
grange reconstruction aR^ M_,M^,e{^) dl lab has M - 1 extrema, where aj^^ ^ ^ ^(^) = 0, one in each interval be- 
tween 2 consecutive roots ^Ri,M-,M+,i,n (I24a| l. Indeed, for any nonzero polynomial p{^) e with M distinct 
real roots we know, by Rolle's theorem lll9[ pp. 215-216], that there is a point where p'{^) = in each of the 
M - 1 intervals between 2 consecutive distinct real roots, these M - 1 points being exactly the M - I roots of 
p'(^) 6 Mm-i[^]. Both aR„M_,M„f(^) 6 Rm[^] (HB, by Proposition |331 and a,,M_,M„((^) & Km[^] by (|23j, 
have M real distinct roots. Therefore, a^^ ^ ^ ^(^) E IRM-i[f] and ^ ^(^) e ]RM-i[f] have M - 1 real distinct 
roots, corresponding to the M - I extrema of aR^ M_,M^,eiO and aiM_,M+AO' respectively. Since (Proposition I3.lt 
aRuM^,M+AO = [R(ij)(ai,M^,M+,()W, by LemmairU o'^^^m^^m^/^) = [^(1;1)K,m_,m+/)](^)' s° correspond- 
ing M - I distinct real roots, which are also the corresponding extrema of aR^ M-,MtAO - [^(1;1)(Q'/,m_,m+/)](^), are 
distant by < ^ (Lemma [l!4l l. □ 

Proposition 3.7 (Factorization of ai^M-,M^AO (I12ab and q;r,,m_,m+/(^) dllab ). Assume M± e Z : M := M_ +M+ > 0, 
and ( £ {— M_,--- ,M+}. Then the fundamental polynomials of Lagrange reconstruction or^^M-M+A^) dllab anc/ 
interpolation ai,M_,MtA^) dl2ab on f/ze stencil Si,M_,M+ (Definition II. lb can i>e factorized as 



o:r,,m_,m+ 



...(f) ^C-D^^'^^j^L^l, ) fl (^-^«..M-,M.x«) (25a) 

"<«='-'>'*s7UV)n«-"' <^''' 

^ ' n~-M_ 

where ^RyM-M*,i,n e M (n e {-M_, • ■ • , M+} \ {^}) are the M real and distinct roots ofaR^^M-,M+AO (Proposition l3.5b . 
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Proof. Every polynomial p(x) e M.[x] can be factorized as p(x) - coeff[;ic''''s*''',p(;ic)] ~ ^pj^ where 

jcp,, 6 C (n e {I,-- - ,deg(p)]) are its deg(p) e N roots (l^, pp. 284-285]. We know that deg(Q'R,_M_,M+/) - 
deg(a/,M_,M,/) = M_ + M+ = M ( fT9cb . The M roots of ff/.M^.M+z© are n e ({-M_, ■ ■ ■ ,M+) \ {/")) c Z and 
the M roots of or^^m-M+AO^ ^Ri,M-,M+,t,n e K. (n e ■ ■ ■ ,M+] \ [I]) are real (Proposition 13. 5t . Furthermore 

□ 



coeff[^ , Q'fi,,M_,M+,f(^)] - coeff[^*',a'/_M_,M+/(^)] are given by ( |19dt . These facts prove (125 

Example 3.8 (Fundamental polynomials Q'/,M_,M+,f(^) ( I12al ) and aRi,M-,M^,t{S) ( 1 11 al l). Consider the fundamental poly- 
nomials of Lagrange reconstruction, as, ,m-,m^AO i ll lal ). and the corresponding fundamental polynomials of Lagrange 
interpolation, aiM_,M+AC> flZaJl, on the stencils (Definition ll.il ) 8,33 (Fig.lTJ and S,,3,4 (Fig.|2]). We know that the 
corresponding polynomials, a'R,,M_,M+,f(^) 6 dl lab and Q'/,M_,M+,f(^) 6 ( I12al ). V/" 6 {-M_, ■ ■ ■ ,M+), have 

M :- M_ + M+ distinct real roots (Proposition l3.5K each root of the fundamental polynomial of Lagrange reconstruc- 
tion ^Ri,M-,M+,e,n ( I24ab being close to the root n € {-M_, ■ ■ ■ , M+] \ {{} of the corresponding fundamental polynomial 
of Lagrange interpolation (I23ab . v/z |^^,,M_,M+/,n -«l < 5 (I24ab . Furthermore (Remark l3.6t both a/,M_,M+,f(^) (I12ab and 
ctRi,M-M+AS) dl lab have M- 1 corresponding extrema, again distant < 5. For these reasons the shapes of ai,M-,MtAO 
(I12ab and Q'r,,m_,m+/(^) (II lab are quite similar (Figs. [T] |2}. For the stencil 8,33 (Fig.[T]i which is symmetric around 
^ = 0, we observe that ^r,33,/',„ i 'LSC e {-3, +3} Vn e {-3, ■ ■ ■ , +3} \ {{}. On the contrary, for the stencil S/3,4 
(Fig. lU which is symmetric around ^ - 5, we observe that there are two integer roots, ^r, 3,4 -3,1 = -i-l e Z and 
fRi,3,4,+4,o = e Z. Although we have not worked out a formal proof concerning integer roots, we can formulate the 
following conjecture (Result l379] l. obtained using symbolic computation. □ 



Result 3.9 (Integer roots of aR,,M_,A/+/(^) dHab ). Let 

M± e {-20, ■ ■ ■ , -H20} : M := M_ + M+ > 1 

Then for M even 

f V^€ {-M_,-- - ,M+} 



2k; ken ^ £,R^M-Ma,n i ^ 



and for M odd 



M = 2k+l; keNo 



V n e (-M_ 



v-rf 1 = -^^^ - if] - 



(26a) 



(26b) 



M+) \ K) 

M+ - M_ -H 1 
2 

M+ - M_ - 1 



(26c) 



^Ri,M-,M^,C,n i ^ 



V f€ {-M_ + I,-- - ,M+ - 1} 

Vn6{-M_,--- ,M+]\[e] 



Verification. By Proposition 13.51 we know that all of the roots of the basis polynomials aR^^M_,M^AO are real. 
Since by ( lllab deg(aRj_M_,M+,c) = M (M :- M_ + M+), there are M real roots, with exactly 1 root in each of the M 
open intervals (n - 5, n -H 5) V n e {-M_, ■ ■ ■ , M+} \ {£] ( I24ab . Hence, if ur^^m^m+A^) has integer roots, these must 
belong to the set {{-M_, ■ • ■ , M+} \ {(]} c Z, /e if Qr^^m-.m+AO has integer roots these must lie on the points of the 
stencil S,,m_.m+ '■- [i - M^, ■ ■ ■ ,i + M+} (Tc\ . except the point / + ( itself. As a consequence, the result was obtained 
by direct calculation of aRi,M-,M+AS) ( II lab 'it, n E {-M_, ■ • ■ , M+] for the range of stencils studied. □ 



Notice that by compaiison of )25bl with 



e+M+ 



1 



M 



M\ \e + M_ 



as can be easily verified by direct computation. 
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5«i 


,3,3,- 


-3,-2 


~ -2.1780 




,3,3,- 


-3,-1 


~ -1.0958 




,3,3,- 


-3, 


~ -0.0307 




,3,3,- 


-3,-H 


~ +1.0307 




,3,3.- 


-3,4-2 


~ +2.0958 




,3,3,- 


-3,4-3 


~ +3.1780 



4 


3 -2 


-1 





12 3 4 






3,3. 


0.-3 


~ -3.1658 








3,3. 


0.-2 


~ -2.0699 








3,3, 


0,-1 


~ -0.9714 








3,3, 


0,-H 


~ +0.9714 








3,3. 


0.+2 


~ +2.0699 








3,3, 


0,+3 


~ +3.1658 





— 1 


\l= o| 






f 





€Bi,3,3,-2,-2 ~ -3.1179 
5fll,3,3,-2,-l ~ -1.1277 
3,3,-2, ~ -0.0446 
«ft,3,3,-2, + l ~ +1.0231 
5fll,3,3,-2,-(-2~ +2.0914 
$fll,3,3,-2,+3 ~ +3.1756 



€fli,3,3,- 


-1,-3 




-3.1533 


?fli,3,3,- 


-1,-2 




-2.0318 


5fli,3,3,- 


-1, 




-0.0796 


5iii,3,3,- 


-1,-H 




+1.0087 


?iSi,3,3,- 


-1,4-2 




+2.0841 


5bi,3,3,- 


-1,4-3 




+3.1719 




■ afli,3,3,f(5) 
a/,3,3,f(?) 

?fli,3,3,f,ii 





,3,3,4-1,-3 




-3.1719 




,3,3,4-1,-2 




-2.0841 




,3,3,4-1,-1 




-1.0087 




,3,3,4-1, 




+0.0796 




,3,3, 4-1, +2 




+2.0318 




,3,3, 4-1, +3 




+3.1533 




fill ,3,3,4-2,-3 




-3.1756 


5fli,3,3,4-2,-2 




-2.0914 


5fll,3,3,4-2,-l 




-1.0231 


,3,3,4-2, 




+0.0446 


5fli,3,3,4-2, + l 




+ 1.1277 


€fli,3,3, + 2,+3 




+3.1179 





I = +2 









,3,3,4-3,-3 




-3.1780 




,3,3,4-3,-2 




-2.0958 




,3,3,4-3,-1 




-1.0307 


5ii 


,3,3,4-3, 




+0.0307 




,3,3,4-3, + l 




+1.0958 


5fl 


,3,3,4-3, +2 




+2.1780 



e = +3 



-4 -3 -2 -1 1 2 3 4 5 -4 -3 -2 -1 1 2 3 4 5 -4 -3 -2 -1 1 2 3 4 5 

Figure 1: Fundamenatal polynomials (Proposition 12. 2t of Lagrange interpolation, ai,M-,M+,e(.^) )12aK and reconstruction, aR^.M-M+.l^^) H lat . 
on the stencil 8,33 (C e {-3, ■■■ ,+3}), and locations of the 6 real roots of each ff^i^^.K?) £ ^6l?] (Proposition 13. 5K ^Ri,M-,M^,e,n 6 
{-M_, ■ ■ ■ , M+l \ {{\), appearing in the factorization (|25a) of a^^ ,M-,M+,f(f) (Proposition l3.7t : notice that in the present case (M := M_ + M+ = 6) 
^Ri,3,3,f,,, « Z Vf e {-3, +31 Vn e 1-3, ■ ■ ■ , +3) \ {H (Result[ID. ' 



3.3. Some identities concerning the fundamental polynomials of Lagrange reconstruction 

To build (SjHi the recursive construction of the weight-functions lUllslIilll o-R^^M-,M+,K^,kS^) for the combination ( fSal l 
of the polynomial reconstructions pRt,M--k^,M+-K,+kS^i + ^Ax; x/, Ax; /) on the substencils (Definition 1 1.2b of S,;m_,m+ 
(Definition ll.il ) to the polynomial reconstruction pR^,M-,M+i^i + ^^x; x,, Ax; /) on the big stencil S,;m_,m+, we will first 
examine (Lemma |4. 2b the elementary subdivision of S,m_,m+ (Definition II. lb into the substencils S,;m_-i,m+ (which 
omits the leftmost point / - M-) and S,,m_.m+-i (which omits the rightmost point ; + M+). To obtain the general result 
for crs,,M_,M+,i,o(^) and o-r^_m_,m+,^,^^^) (Lemma |4!2] | we need to show that the leading terms of the approximation 
error (l8ct of the Lagrange reconstructing polynomial, on 2 overlapping stencils of equal length, but shifted by 1 cell. 
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.3,4,- 


-3,-2 




-2.1897 


Sfll,3,4,- 


-3,-1 




-1.1131 




,3,4,- 


-3, 




-0.0538 




,3,4.- 


-3,-H 




+ 1 




,3,4,- 


-3,4-2 




+2.0538 




.3.4,- 


-3,4-3 




+3.1131 




.3.4.- 


-3,4-4 




+4.1897 



e = -3 



5fli,3.4, 


-2,-2 




-3.1289 


5iJi,3,4, 


-2,-1 




-1.1437 


5fli,3,4, 


-2, 




-0.0674 


5fli,3,4, 


-2,4-1 




+0.9925 


5fli,3,4, 


-2, +2 




+2.0492 


?fli,3,4. 


-2, + 3 




+3.1102 


€-Rl,3,4, 


-2, +4 




+4.1881 



€iii,3,4, 


-1,-3 




-3.1639 


5hi,3,4, 


-1,-2 




-2.0464 


5-Ri,3,4, 


-1, 




-0.1012 


5fli,3,4, 


-1,4-1 




+0.9782 


5-Ri,3,4, 


-1,4-2 




+2.0417 


5fli,3,4, 


-1,4-3 




+3.1059 


(Ri.3A, 


-1,4-4 




+4.1857 




5fll.3,4, 


0,-3 




-3.1761 


5fll,3,4, 


0,-2 




-2.0846 


Sfll,3,4, 


0,-1 




-0.9892 


5fll,3,4, 


0,4-1 




+0.9419 


5iil,3,4, 


0,+2 




+2.0271 


5fll,3,4, 


0,4-3 




+3.0987 


Cfil,3,4, 


0,+4 




+4.1821 



3,4,4-1, -3 




-3.1821 


5fll,3,4,4-l,-2 




-2.0987 


3,4,4-1,-1 




-1.0271 


3,4,4-1, 




+0.0581 


3,4,4-1, +2 




+1.9892 


5iJl,3,4,4-l,+3 




+3.0846 


?-R.l,3,4, + l,+4 




+4.1761 




■ afli,3,4,f(5) 

a;,3,4,<'(0 

?fll,3,4,f,ii 



4 


3-2-1 1 


2 3 4 


5 -4 


3-2-1 


1 


2 3 4 




5i!i,3,4,4-2,-3 ^ 


-3.1857 




3,4,4-3, -3 




-3.1881 






Ciii,3,4,4-2,-2 ^ 


-2.1059 




3,4,4-3, -2 




-2.1102 






3,4, 4-2,-1 ~ 


-1.0417 




3,4,4-3,-1 




-1.0492 






5iJl,3,4,4-2, ^ 


+0.0218 




^fli, 3,4,4-3, 




+0.0075 






5fll,3,4,4-2, + l ^ 


+1.1012 




3,4,4-3, + 1 




+ 1.0674 






5fll,3,4,4-2,+3 ~ 


+3.0464 




€iJi, 3,4,4-3, +2 




+2.1437 






5Hi,3,4,4-2,+4 ~ 


+4.1639 




3,4,4-3, +4 




+4.1289 






5hi,3,4,4-4,-3 




-3.1897 


5-Ri,3,4,4-4,-2 




-2.1131 


Sill, 3,4,4-4,-1 




-1.0538 


5-Rl, 3,4,4-4, 







5fll,3,4,4-4,4-l 




+1.0538 


5Hi,3,4,4-4,+2 




+2.1131 


5fli,3,4,4-4,+3 




+3.1897 



1 


t = +i 




/ 


£. — - 



-3-2-10123 

Figure 2: Fundamenatal polynomials (Proposition 12. 2) of Lagrange interpolation, a/.M-,Af+.f(f) )12aK and reconstruction, a«i,jM-,M+,f(f) Ulat , 
on the stencil 8,34 {t e j-3, ■■• ,+4|), and locations of the 7 real roots of each Q'R|.3,4/(f) e K7[f] (Proposition 13. 5K ^R^,M-Mt,t,n (" £ 
{-M_, • ■ ■ , M+] \ {/)), appearing in the factorization (|25a} of alt^ ,M_,M+,f(?) (Proposition l3.7t : notice that in the present case {M := M_ + M+ = 7) 
^Ri, 3,4,-3,1 = +1 e Z and ^R,, 3,4,4-4,0 = e Z (Result|32). 



are different. Since the error-expansion (l8ci polynomials ( lllcl i. /lfi,,M_-i,M+,M(^) and /ifi,,M-,M+-i,M(^), are of degree 
fn8, p. 294, Proposition 4.7], they can be projected on the basis {q;R],m_,m+,^(^)> ^ £ {~M^, ■ ■ ■ , M-i-}) of RM[f] 
(Proposition l3.3b . and the same projection is possible for the polynomials {Q's,,M_-i,M+,f(^), ^ £ {-M- + 1, • • ■ , M+}} € 
Rm-i[^] c Rm[^] and {aR„M_,M,-i,f(^), ^ e {-M-, ■ • ■ ,M+ - 1)) e Rm-i[^] c Rm[^]. ' 

Proposition 3.10 (Identities on Q'R,,M_,M+,f(^) (II lal l and /iR,,M_,M+,«(^) (II lal l). Assume the conditions and definitions 
of Proposition 12.21 and consider the stencil S,,m_,m+ (Definition II. lb and its substencils (Definition 11.21 ) S,,m_,m+-i 
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and Si 



4^. The following identities hold E 



Q'fii,M--l,M+,M+(^) -(-1) (XRuM-M^-I-mSO 

aR^M-M,-\M) *aRuM^-\,M,M) ^ ^ 6 {-M_ + 1, ■ ■ ■ , M+ - 1} 



A. 



Ri,M-M+-l,M 



fT9bl 



(27a) 
(27b) 
(27c) 
(27d) 

(27e) 



Proof. Let p{x) e ]Rm[-^] be a polynomial of degree < M. Then by |0, Theorem 5. 1, p. 296] its reconstruction pair 
q{x) := [^(i; Av)(p) ](x) 6 Km[-^], and deg(^) = deg(p) (Remark|23|i. 

Proof of ( I27al l: By Proposition 12.21 V/7(x) € Mm-i[-^], its reconstructing polynomials (Definition 11.4) on the 2 
stencils S,_m_,m+-i '■- {i - M^, ■ ■ ■ ,i + M+ - 1) and S;,m__i_m+ := {i - + 1, ■ ■ ■ , ; + M+), which contain the same 
number of M points but are shifted by 1 cell, are exactly equal to the reconstruction pair of p(x) (Definition ll.4l i 
q(x) :- [R(\;Ax)(p)](x), because of dSbl l. since p(x) e Rm-iM => p'^^Hx) - Or^_j[x](x) > M. Hence, by (l8at . 



M+-1 



= P«i,M_-i,M+(JC; + ^Ax; x,-. Ax; /?) = 



'i*q'(x; + ^Ax) [7?(I;Ax)(p)](^; + ^Ax) 



P«i,M_,M+-i(-Xi + ^Ax; Xj, Ax; p) ^ aKi,M_,M+-i/(^) p(-^i + ^ Ax) 

cslM^-lm+.K^) /5(jc/ + £ Ax) 

f=-M_ + l 

Vp(x)€IRm-iW 

Vx,- 6 R 
VAx e R>o 



whence 



a«,,M_,M+-i,-M_(^) p(xi - M_ Ax)+ 
^ (aiJ,,M_,M+-i/(^) - a«i,M_-i,M+,f(^)) /'(■^i + f Ax)- 

Q'«,,M_-l,M+,M+(^) P(-^,- + M+ Ax) = 

Applying (I28bl i to the polynomial 



V;?(x)eRM-iW 
V^eR 
Vx; e R 
VAx e R>o 



I [-^l 3 Y\ ^' ~ - Vx G [xi - (M_ - l)Ax, ■ ■ ■ , X; + (M+ - l)Ax} 



m=-M_+l 



yields 



q;«,,m_,m+-i,-m_(^) 



]~[ (-M - m) 



Vm=-M_ + 1 



Vm=-M_ + 1 



ff«i,M_-l,M+,+M+(^) 



=0 V te(-M_ + l,--,M+-l) 

Yl i+M^-m) 

Vm=-M_ + 1 



V^e R 
VAx e R>o 



(28a) 



(28b) 



(28c) 



(28d) 
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le 



OlRx,M-.M+-l-M-{^) 



\ ( Mt- l 

- Q'R,,M_-l,M+,+M+(^) 



Vm=-M_ + 1 



P] (+M+-m) 



^m=-M_4-l 



e R (28e) 



Since 



M+-1 



Y[ (+M,-m)'"=^-'" f[ /t= (M-1)! 

m=-M_ + l it=M-l 
M+-1 -(M-1) 

PI ^_M^ W /t (M-1)! 



m=-M- + l 



k=-l 



using (l28fl[28g]l in (l28d] i proves ( l27al i. 

Proof of ( I27bb : Applying ( I28bl ). successively for A: e {-M_ + 1, ■ ■ • , M+ - 1}, to the polynomials 



m=-M_ + l 



yields 



m=-M- + l 

\fke {-M_ + I,-- - ,M+ - 1) 



m=-M- + l 



and using (28g 



(-l)^M! UM 

a«,,M_,M+-l,i(^) - Q'«,,M_-l,M+,/:(^) a'R|,M_,M+-l,-M_(^) ORM-l[f](^) 



n - 



y ke {-M_ + I,-- - ,M+ - 1) 



proving (l27bli bv (fT9bliPI 



Proof of (|27cl[27dll: Notice first that by (fTTcT l. for n = M + 1 = M_ + M+ + 1, 



f=0 



,M_,M+,M+l(^)— - 5 - ^ - 5) 



-A'iJi.M-.A/t.M+lC^) 



VM+ e Z : M := M_ +M+ > 1 



(28f) 
(28g) 



3 ]~~[ (jc - Xi - mAx) -0 Vx e {x; - (M_ - l)A;ic, ■■■ ,Xi + M+Ax] \ {xi + kAx] (29a) 



(29b) 



(29c) 



(30a) 



in=-M- + 1 »=^+M_ - 1 
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where fiR,,M_,M+,n(^) is defined by ( lllbl i. By ( I30al ) 

^Ri,M_,M+-l,M(,gJ - A'«,,M_,M+-l,M(,gj 
J ^^^^^ ^er^ / ^ 

Since \/p(x) e M.m[x] ^ p^"\x) = Or„[.i.](x) \ln>M+ l,we have by (l8al[8b]i. 



± eZ: M:=M_+M+ >2 



aRx,M-,Mt,e\^ P(xi + 1 Ax) = ^(x,- + ^Ax) 



{ ^p(x) e Rm[x] 

qix) = [^(l;Ax)(;')]W (O 

VAx € M>o 



^ aRuM^-\,M.M) P(xi + f Ajc) ^^q{xi + ^Ax) + IJ.r,M--\,m.,m(^)^x^ P^'^\xd 



=-M_-l 



for the reconstructing polynomials Pr,,m_,m+-i(jC( + ^Ax;xi,Ax;p) (I30dl l. pR,,M_,M+ixi + ^Ax;xi,Ax;p) 
PRi,M--i,M+(xi + ^Ax; Xi, Ax; p) dSOfl ). Consider the polynomials 



M+-1 



■m[-^] 9 ]~~[ (-^ - -JC/ - otAx) =0 Vx e {x; - M_Ax, ■ ■ ■ , x; + (M+ - l)Ax) 

m=-M_ 

[x] 3 (x - x; - otAx) =0 Vx e {x,- - (M_ - l)Ax, ■ ■ ■ , x,- + M+Ax} 



Obviously, 



dx^ 



(x - Xj — mAx) 



\fl;=-M_ 



(x — X, - mAx) 



M! \ Vxi e R 
I VAx6K>o 



(30b) 
(30c) 



(30d) 

(30e) 

(30f) 
and 

(30g) 
(30h) 



(30i) 



Applying ( I30eb to the polynomials ( 30g I30hb . (I30db to the polynomial ( 30g i, and (I30fb to the polynomial (I30hb . we 
have, using (I30ib . 



M+-1 



([30e]|30g| ^ aR,,M_,M,,+MM) A^*' ]~[ (M+ - m) = R(i,^^) fl ~ ^' ~ "^^^^ 

m--M- \m--M- 

(130eJ|30a ^ a«„M_,M„-M_(a Ax^ ]^ (-M_ - m) = 7;(1;a,) (x - x, - mAx) 



(130d||30gb^0^j;(i;A.,) 



(|30B|30hli ^ = 



(x — X, - mAx) 



Vm=-M_ + 1 



+ ;U«,,m_,m+-i,m(^)Ax'^ M! 



(I;Aa:) ]~~[ (x-X,-mAx) 

VxeM ^Ax:=x-x,- Vx/ e R VAx e R>o 



+ ;"«,, m_-i,m+,m(^)Ax*^ M! 



(30j) 
(30k) 
(301) 
(30m) 
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and combining pOj[) with ( 13011 ), and ( I30kl ) with ( 130ml ). we have 



aR,,M_,M„+M,© n ~ '"^ ~ J"«,,M_,M,-1.M© M\ E R (30n) 

m=-M_ 



,m-,m.,-mA^) n (-^--'")'^^-J"«.,M-i,M„M©^! V^eM (30o) 



whicnUby (l30b][30c]i prove dZTcllZTdl l. 

Proof of ( I27eb : Applying (I30fb to the polynomial pOg[) yields 



Q'«i,M_i,M^.,+M^.(^) Ax*^ ]~~[ (M+ -m) '^=^i;(i;Ax) ]~~[ (jc - X/ - otAx) 

m=-M_ 

Combining ( l3Tal ) and ([Ml) yield^El 



+ //r,,m_-i,m,,m(^)Ax^M! (31a) 



\m=-M- 

VxeK ^Ax:=x-x; Vx;€R VAxeMx) 



ffSLM-LM+.+M+C^) '"^^ -/^iJi,M_,M+-l,M(^) +/^iJi,M_-l,M+,A/(^) £ M (31b) 

which by ( l30bl [30cT l proves (l27el l . □ 

4. Reconstruction by combination of substencils 

4.1. Substencils of Si^M-,M* 

WENO reconstruction on S,;m_,m+ achieves high-order in smooth regions and monotonicity near discontinuities 
by a nonlinear (depending on the values fi+i of fix) on the points of the stencil S,-^f m+) combination of recon- 
structions on substencils whose union equals the stencilOU Central to this development is the determination of the 
underlying optimal (linear in / in the sense that the weight-functions depend only on x and not on /) combination of 
the reconstructing polynomials on the substencils to exactly obtain the reconstructing polynomial of the entire stencil. 



Example 4.1 (Substencils (Definition ! 1.21 )). Notice that a given stencil S, m_,m+ can be divided into different families 
of substencils (Definition 11.21 ). depending on the chosen value of Ks < M - \ ( l2b] i. The 0-level of subdivision 
{K^ - 0) corresponds to the original stencil, without subdivision. The (M - l)-level of subdivision {K^ - M - Y) 
corresponds to the subdivision of the original stencil to A's -H 1 - M substencils of length equal to 1 cell, ie to the 
substencils [S /,m_.m_+i , • ■ • ,5' ,.m+-i,m+ ), on each of which polynomial interpolation, and as a consequence polynomial 
reconstruction (Remark I2.4K are of degree 1 (linear). As an example, we consider the successive subdivisions of the 
stencil S, 3 3 (Fig. O which corresponds to a stencil symmetric around point /, and of the stencil 8, 3 4 (Fig. |4|i which 
corresponds to a stencil symmetric around point i + \. We called the substencils of Definition [L2] Neville substencils 
lfl2ll because they are those used in the Neville algorithm pp. 207-208] for the recursive construction of the 
interpolating polynomial. □ 



By analogy with )28fl|28g) we have 

M+-1 1 M+ -M 

f] (+M,-m)* = =*-'"f]/t = M! ; f] (-M_ - m) f] yt = (-1)" M! 



m=-M- k=M m=-M- + 1 

Shu f3,2l uses the tenns big stencil and small stencils to denote the stencil and its substencils. 
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if. = (il/ - if, + 1 = 7 points, M -1 if, = 6 intervals) 



if, = 1 (J\i - if, + 1 = 6 poiiits, M ■ 



if, = 2 {J\i - A", + 1 = 5 points, il - 



A", = 3 (Jli - if, + 1 = 4 points, ,1i -! 



A, = 4 [M ~ if H + 1 = :j points, M -i A', = 2 intervals) 

• • • ■ 



A', = 5 [M - if, + 1 = 2 points, A[ -, if, = 1 interval) 



-e- 



A, = 5 intervals) 



if, = 4 intervals) 

• e o 



if, = 3 intervals) 







Si, 2, 2 



S.,1,2 
^1,0, 3 



Si,r,r 

Si,0,2 
S,,-1.3 



ure 3: Successive subdivisions of the stencil S, 3 3, for different values of A'j e |0, ■ ■ ■ , M - 1 = 5) (Definition ! 1.2t . 



A', = (M - A', + 1 = 8 points, M - K\ = 7 intervals) 
A, = 1 (i\i - if, + 1 = 7 points, 1\1 - A', = 6 intervals) 



A, = 2 {M - A, + 1 = 6 points, M - A, = 5 intervals) 
• • • • • • 



A, = 3 (Ai - A, + 1 = 5 points, M - A, = 4 intervals) 

• • • • • e e- 



A, = 4 (i\i - if, + 1 = 4 points, ili - A, = 3 intervals) 

• • • • — < — e e e o ^.ao 

• • • • e e o s,.2,i 



if, = 5 {M - if, + 1 = 3 points, M - A, = 2 intervals) 

• • • e — I — e e 



if, = 6 (A/ - if, + 1 = 2 points, M - A, = 1 interval) 

• • e e e e 



s.,o,r 
s.,-r,2 

Si,-2,3 
Si,-3,4 



ure 4: Successive subdivisions of the stencil S, 3 4, for different values of e |0, 
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,M- 1 



= 6) (DefinitionO- 



4.2. (Ks — lyievel subdivision 

The starting point for developing a recursive formulation for the weight-functions is to consider the {K^ - l)-level 
subdivision of S, m_,m+ The resulting substencils S,,m_,a/+-i and S, « have equal lengths of M - 1 cells (M points), 
but are shifted by 1 cell (Figs. |3] |4]i. If a {K^ - l)-level subdivision rule can be established, then it can be readily 
extended to {K^ > l)-levels using the general recurrence relation proven in [ 12, Lemma 2. 1]. 

Lemma 4.2 (Rational weight-functions for {K^ - 1 )-level subdivision). Assume the conditions and definitions of Propo- 
sition l2!2l and consider the stencil S,;m_,a/+ o.nd its substencils (Definition 1 1.2b S,;m ^f^-i and S,yv/ Define the 

functions crj?,_M_,M+,i,o(^) and crs|,M_,M^_,i,i(^) by 

Q'«i,M_,M+,-M_(^) {27] ^R,,M_-I,M^_,m(^) , 

o-R,M-,M,xm> '— TTT - 1 77\ ] TTT ^^^^> 

.^s _ fffii, M_,M+,M+(^) Ig} '1j?,,M-,M+-1,m(^) /tiu\ 

o"Ri,M_,M+,i,i(?) :- — - — — (J^o) 

satisfying the consistency condition 

CR,,M_,M+,1,o(^) + 0"R,,M_,M+,l,l(^) = 1 (32c) 



Then the reconstructing polynomial on Si^M_,M+ (Proposition I2.2l l can be constructed by combination of the recon- 
structing polynomials on the 2 (Ks — V)-level-subdivision substencils as 



(32d) 



PRi,m_,m^ (xi + f Ax; Xi, Ax; f) = 

_ AR^,M--\M^,MiS) PR,.M-,M^-l(Xi + ^Ax; Xj, Ax; f) - /iRi,M_,M^-l.M(^) PRx,M--\,MS^i + ^^^'^ ^i' ^^'^ /) 

'^Ri,M_-1,M+,m(^) - ^Ri,M_,M+-1,m(0 

E M Vx; E M VAx £ Kx) V/ : M R 

and can be represented, almost everywhere, as 

PRt,M_,Mjxi + ^Ax; Xi,Ax;f) ^o-r^^m^^m^^iM) PR^M^.M^-iixi + ^Ax;xi,Ax;f) 

+o-Ri,M-,M*,iM) PRi,M--\,MSxi + ^Ax; Xi, Ax; f) (32e) 

V^E R\ R : ffR,,M_-i,M„M,(^) = 0} Vjc;eR VAxeMx) V/ : R — » M 
The functions o-r^^m_^m^,i,o(^) ond cr^^ m_,m+,i,i(0 satisfying ( I32cll32e] l are unique. 
Proof. By Proposition [TTO] we have 

Q'«i,A/_,M+,-M_(^) l27all27cl (-I)*'"' '^R,,M_-1,M+,m(^) {27e) 't^i,M_-l,M+,M(^) 



«fii,M_,M+-l,-M_(^) (-1)''^ ' Q'/fi,M_-l,M+,+M+(f) ^R,,M_-\,M+,MiO - ■Ir,,M_,M^-\,m(0 

proving (I32al l. and 

(XRiM-M^mSO I27dlj27el -^Ri,M_,M+-1,m(^ 
Q'«,,M_-l,M+,M+(^) ^Ri,M_-I,M+,m(^) - 'ifi,,M_,M+-l,M(^) 



(33a) 



(33b) 



proving (I32bl l. Obviously ( I32cl l holds because 



'^Ri,M_-I,M+,m(^) - '^Ri,M_,M+-1,m(^) 

-'iR,,M_,M+-l,M(^) 
^Ki,M_-1,M+,m(^) - ^Ri,M_,M+-1,m(^) 
_'lfii,M_-l,M+,M(^) - ^Ri,M_,M+-1,m(^) 
'^Ri,M_-I,M+,M(f) - 'iRi,M_,M+-l,M(^) 
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= 1 (33c) 



Proof of (l32e]i: Let p(x) e Rm[x] (O. Then, by (UH), q(x) [^(1;Ax)(/?)](.x) e Mm[.x]- Since Vq(x) e Mm[.s:] 
^'"'(jc) = Or„[x] Vn > M + 1, we have, by application of (l8cb . and taking into account Remark lZ4l 



PR^,M_,M^-i(xi + ^Ax; Xi, Ax; p) '^^q(xi + ^Ax) + /tj?,,M_,M+-i,M(^) Ax'^ + ^Ax) 



(34a) 



PRt,M-,MMi + ^^x; Xi, Ax; p) ^^q{xi + ^Ax) 



q(x) := [R(UAx)(p)](x) (O 
V^€ R 
V;c,- e R 
VAx € M>() 

PR,,M_-i,M,te + ^Ax; X/, Ax; p) ^^q(xi + ^Ax) + /Ir,,m_-i,m+,m(^) Ax*' ^^'^(x/ + ^Ax) 
Combining (I34al i weighted by ( I32al i.and ( I34cl i weighted by ( I32bb yields 

o"R,,M_,M+,i,o(^) ;'R,,M_,M+-i(-Xi + ^Ax; Xi, Ax; p) + crR,_M_,M+,i,i(^) P«i,M_-i,M+(-Xi + ^Ax; Xi, Ax; p) = 



(34b) 



(34c) 



(o"R,,M_,M+,l,o(^) + crR,,M_,M+,I,l(^)) ^(-^i + ^Ax) 
+ (o"R,,M_,M+,l,o(^) /1«,,M_,M+-1,m(^) + CrR,,M_,M+,l,l(^) '^R,,M_-1,M+,m(^)) Ax*' g'<'^'(x,- + ^Ax) 



"-^-•^(x,- + ^Ax) 



/1«i,M_-1,M+,m(^) ^R,,M_,M+-l,M(f) -^R,,M_,M+-1,m(^) '1ri,M_-1,M+,m(^) \ . fMW >- » x 

Ax q '{Xi + ^Ax) 



'lfii,M_-l,M+,M(^) - 'tRj_M_,M+-l,M(^) '^Ri,M_-1,M+,m(^) " 'l/fi,M_,M+-l,M(^) 

f Vp(x) e RmIx] 



t34bl 

=^(xi + ^Ax) := PriM-mMi + ^Ax; x,-. Ax; p) 



q{x):^[R(,,L,){p)]{x)m 
V^eR 
Vx; e R 
VAx e R>o 



(34d) 



^cr«j3^ _Af+,i,o(^) P«,,M_,M^-i(-Xi + ^Ax; Xi, Ax; p)+ 

o'«,,M_,M+,i,i(^) P«,,M_-i,M^(-«i + ^Ax; Xi, Ax; p) = pR^,M_,MMi + ^Ax; x;. Ax; p) (34e) 
V/7(x)€Rm[x] Vf€M\{^eM:cKR,,M_-i,M,,M,(^) = 0) Vx,- 6 M VAx 6 M>o 

which shows that (I32el l is vaUd V/(x) e Rm[-«]- Using the representation (ISal l of the reconstructing polynomial in 
(134^ 

'^'^(q;«i,M_,M+,-M_(^) - O-R,_M_,M+,l,0(f) Q'«i,M_,M+-l,-M_(^)) p(JCi - M_ Ax) 

M+-1 

+ ^ (Q'Ri,M_,M+,f(^) - 0"«],M_,M+, l,o(^) 0:Ri,M-,M^-iA^) ~ 0"«i,M_,M+,l,l(^) Q'R,,M_-l,M+,f(^)) /'(-^i + t Ax) 

f=-M_ + I 

+ (a«,,M_,M+,+M+(^) - 0"R,,M_,M+,l,l(^) Q'«i,M_-l,M+,+M+(^)) P(jC/ + Ax) 

mm 

= 2_i {^R\M-M^A^) ~ O"«i,M_,M+,l,0(^) Q'Ri,M_,M+-l,f(^) - 0"R,,M_,M+,l,l(^) Q'R,,M_-l,M+,f(^)) /'(-*^! + ( Ax) 

Vp(x)eRM[x] e R \ e R : aR„M_-i,M„M,(a = 0} Vx; e R VAx € R>o (34f) 



where we used 



Q'R,,M_,M+,-M_(^) '^CrR,,M_,M+,I,o(^) CR, ,M_,M+-1,-M_ (^) 6 
<3'Ri,M-,M+,M+(^) '^CrRi,M_,M+,I,l(^) <3'Ri,M--I,M+,M+(^) V^ € 1 



(34g) 
(34h) 
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Applying ( I34fl ) successively to the polynomials 



dx] B Y\ (x-Xi-mAx)-0 Vx e {x,- - M_Ax, ■ ■ ■ , jc,- + M+Ax} \ {x/ + Mx} (34i) 

yields 



m=-M- 



0!Ri,M_,M^,e(^) -O"fi,,M_,M+,l,0(^) <3'fi,,M_,M+-I,f (^) 

Combining the representation (l8at of the reconstructing polynomial with p4g|l34h1|34j| i proves ( I32eb . V/ : M — > K. 
Proof of ( l32db : Obviously, the functions cr^,,M_,M+,i,o(^) and crgj_M_,M+, i,i(^) are defined everywhere (Vf e M) 



except at 

{&,,M_-i,M+,M+,n; n 6 {-M_ + 1, ■ ■ ■ ,M+ - 1} (35a) 

l24l 

W 6 M : aR,,M_-i,M,,+M,(^) = 0) (35b) 

^{^eR: a«„M_,M,-i,-M_(^) = 0} (35c) 

e M : /lfi|,M_-i,M+,M(^) - ^Ri,M-,M+-i,M(f) = 0} (35d) 

Recall that (Proposition 13. 5l l all of the M - I roots of the polynomial Q'r,,m_-i,m+,+m+(^) are real (|24] |. However, using 
( I27all27bl i in (l32eb vields 

. ^ . .^ (32i][27i]|27b) 

PR,,M-,Mt \Xi + ^Ax; X;, Ax; /) 
_ ^j?i,M_-i,M^,M(g) pR^,M_,M^-\{xi+^Ax; x/,Ax;/) - /lfii,M_,M^-i,M(^) /?j;i,M_-i.Mt(^/ + ^Ax;x,-, Ax;/) ^^^^^ 

'tj?i,M_-l,M+,M(^) - 'lRi,M_.,M+-l,M(f) 

E R \ 6 R : Q'R,,M_-i,M„M,(^) = 0) Vx; 6 M VAx € R>o V/ : M ^ M 



or equivalently, using the representation (l8a> of the reconstructing polynomial in (I35el i 

/'R, ,M_,M+ (Xi + f Ax; X,-, Ax; /) = 

^Ri,M_-l,M+,M(f) aRi,M_,M+-l,-M_(^) ^, , x 

/(X; - M_Ax) 



) 0!R,,M-,M+- 

^Ri,m_-i,m+,m(^) - 'Ir,,m_,m+-i,m{0 



Z'iRi,M--I,M+,M(^) Q'Ri,M_,M+-1/(^) - 'tRi,M_,M+-l,M(^) ff«i,M_-l,M+/(^) „, „. , 
— — /(X; + tAx) 
/IP. fl/f _1 H^il") — /IP. fl/f fl/f _l fl-fl^l 

f=-M_ + l 



+ 1 ] — /(x, + M+Ax) (35f) 

^Ri,M--l,M^,M\^) - 'IRi.M-,M+-1,M(^) 

V^€R\{^€R:aR„M__i,M„M,(^) = 0) Vx/ e M VAx e R>o V/ : R ^ R 
To prove ( I32db by (I35fb we need to show that it is valid e M. Rewriting ( l33all33bT i we have 

(33a)/lR,,M_-l,M+,M(f) Q'Ri,M_,M+-l,-M_(^) . ™ 

fffii, M_,M+,-M_(?) - — — : — e K ^ 

^RuM--l.M+.My?) ~ 'iRi,M-,M+-l,M\?) 

(^^RuM_-\,M+Mi^) ~ ^Ri.M-,M+-1,m(0) I (^Ri,M_-1,M+,m(^) a'Si,M_,M^_-l,-M_(^)) (35g) 

(33a)-/lR,,M_,M+-l,M(f) aR,,M_-l,M+,M+(^) „ 

Q'Ri,M_,M+,+M+(f) 5 jzr V^eK^ 

[^Ri,M--1,M+,m(^) - '^Ri,M-,M+-1.m(^)) I ('^Ri,M-,M+-l,M(f) a;Ri,M--l,M+,M+ (f)) (35h) 
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Recall that by ([TTcl)deg(.lfi,,M_,M,-i,M) = deg(.ts,,M_-i,M„M) = M Proposition 4.7, p. 294], by dTTHl deg(aR,,M_,M+-i,- 
deg(Q;R,,M_-i,M+,Mj = M - 1 and deg(Qffi,,M_,M„-M_) = deg( aR,,M_,M+,+Mj = M, and by (l27e] i deg(^R,,M_-i,M+,M - 

^Ri,M_,M+-1,m) = deg(Q;j;,_M _-l, M+,M +) = M - 1. 

Using again ( 127 al l27bT i in ( 34j i yields 



€ R \ € R : a«,,M_-i,M,,M,(^) = 0) e {-M_ + 1, ■ ■ • ,M+ - 1) 

|27i][22b][34]) /Ix,,M_-1,M+,m(^) <3'R,,M_,M+-1/(^) - ^Rj,M_,M+-1,m(^) Or m_-1,M+A0 .oc-x 

ffSLM^.M^K^) = ' ] tt: — TTT (35i) 

Since the set e R : Q'R|,m_-i,m+,m+(^) - 0) contains only M - 1 isolated points (Proposition [33]l, the result of the 
polynomial division ( I35ib must be valid e M, implying 

V^€ {-M_ + I,--- ,M+ - 1} (35j) 



By ( [35gl[35hl[35j] l, we have that ( |35B is vaHd e M, proving ( l32dl i. 

Proof of uniqueness: We have proved existence by construction, E M \ e R : Q'r,,m_-i,m+,m+(^) = 0), 



of rational weighting functions, crs,,M_,M+,i,o(^) ( I32al i and o"j}|^ ^ i i(^) ( |32bb . satisfying the consistency relation 
( I32cb . which combine the reconstructing polynomials on the substencils (Definition II .2b S,;m-,m+-i and S,;m_-i,m+ 
into the reconstructing polynomial on S, m-,m+ ( I32el) . To prove uniqueness, recall that by (I32cb cr^,,M-,M+,i,i(^) = 
1 - cr«|,M_,M+,i,o(^), and rewrite (I32et as 

PR,,M_,MAxi + X;, Ax; /) t^^lplil a-fi,,M_,M+,i,o(^)(/?R,,M_,M+-i(-^/ + ^Ax; x,-, Ax; /) 

-PRuM^-i.mA^' +^Ax; x;,Ax;/)) 
+/'R,,M_-i,M+(-«i +^Ax;x;,Ax;/) (36a) 

E M \ 6 R : Q'r,,m_-i,m,.m,(^) = 0) Vx,- e R VAx e R>,) V/ : R ^ R 

Hence, assuming the existence of 2 different weight-functions [ctrj^/ j o] (^) [cr^jj^ ^ j oj (^) satisfying ( I36ab 
e R V/ : R — > R would imply /5r,,m_,m^-i(-X/+^Ax; x,-, Ax;/) - /?r,,m_-i,m+(^/+^Ax; x,. Ax; /) e R V/ : R — > 
R. This is obviously a contradiction, since, by Proposition 12. 21 the 2 polynomials Pr,,m_,m+-i(-^! + ^Ax; x,, Ax; /) and 
Pr,,a/_-i,m+(jc/ + ^Ax;x,-, Ax; /) are defined by 2 different sets of values, {/-«_,■• ■ ,/+m+-i} and {/-m_+i,- ■ ■ ,/+mJ, 
respectively^ □ 

Remark 4.3 (( I32db vs ( I32eb ). The expression (l32dl l of the Lagrange reconstructingpolvnomial pg, m m, (x, +f Ax; x, . Ax 
on S i^M_M+ is valid £ R, because the rational expression (I32db yields exactly Pr^,m_,m^ (xi + ^Ax; x,-, Ax; /) by poly- 
nomial division. On the other hand, the weight functions crs| M_,M+, i,o(^) (I32ab and crs,,M_,M+,i,i(^) (I32bl) . are not 
defined at the poles of the rational expressions (I32al I32bb . where the representation (I32eb is not possible. □ 



Corollary 4.4 (Identities for (/Ts = l)-level subdivision). Assume the conditions and definitions of Lemma l4.2l Then 
the following identities hold 

<^R,,M-,Mt-M-(0 -CRuM-.M+.lfliO O!Ri,M-,M+-l-M-(0 (37a) 
as,,M_,M+/(^) =CrS| ^f_,M+,l,o(^) ffRi,M_,M+-l/(^) 

+0'R, ,M_,M+, l,l(^) »Ri,M_-l,M+/ (^) V/'e {-M_ -I- I,--- ,M+ - 1} (37b) 

a'Ri,M_,M+,+M+(^) -0'R,,M-,M+,l,\(^) Q'Ri,M_-l,M+,+M+(^) (37c) 

=Crsi,M-,M+,l,o(^) ^Ri,M-,M+-1,m(^) + CrRi,M_,M+,l,l(^) 'iRi,M_-l,M+,M(^) (37d) 

^R,,M_,M+,n(^) =CrR,,M_,M+,l,o(^) '^R,,M_,M+-l,n(^) + 0"R, ,M_,M+, 1,1 (^) 'iR,,M_-l,M+,n(^) Vn > M -(- 1 (37e) 

E R \ E R : aR„M_-i,M„M,(^) = 0) 



'A more detailed proof is given in Proposition ITT] 
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Proof. We have already proved (I37al i as (34g , (I37cl i as (I34hl i. and (I37bb as (34j i. They are summarized separately 
here for future use. Identity (I37d) follows directly from the definitions of cr^i,M_,M+,i,o(^) ( I32al ) and crs,,M_,M+,i,i(^) 
(I32bl ). and was used in the calculations leading to (I34dl ). To prove the relation (I37e) . we replace pR,,M-,M+(^i + 
^Ax;xi,Ax;f), pRt,M-,M+-i(^i + ^Ax;x;,Ax;/) and pR,,M--\,M+(^i + ^Ax; Xj, Ax; f) in (I32el ) by their expansions in 
terms of the derivatives h^"\xi + ^Ax) (IScT l. and obtain, using (l32cll37dT l 

2 ('^R„M_,M„i,o(^) ^R„M_,M,-i,„(^) + tr«„M_,M„i,i(^) i«„M_-i,M„„(^)) Ax" h^"\x + ^Ax) ™^ 0(Ax^-+') (38) 

n=M+l 

6 K \ 6 M : a«„M_-i,M„M,(^) = 0) \/he C'^-+'(R). 

Using polynomials ^(jf) € R„[x] (Remark |2.4| |. recursively for n > M + 1, in (l38T l proves ( I37eb . by induction. □ 

4.5. (/Ts > lyievel subdivision 

We have shown in lll2i Lemma 2. 1] that if a general family of functions Pm_,mS^)^ depending on 2 integer indices 
M± e Z : M_ + M+ > 1, is equipped with a (/Tj = l)-level subdivision rule, defined by a relation of the form dSal l with 
A's = 1, then, by recurrence, we can construct weight-functions satisfying dSal l V/Tj < M - 1. By Lemma l4~2l we can 
always define uniquely the optimal weight-functions ( |32] | of the {K^ - l)-level subdivision of S, m_,m+ (Definition ! 1.2ll . 
Therefore, the Lagrange reconstructing polynomials Pr^m-mA^i + ^Ax; x,, Ax; /) (ISal i satisfy the conditions of lll2[ 
Lemma 2.1]. 

Proposition 4.5 (Recursive generation of weight-functions for the Lagrange reconstructing polynomial). Assume the 
conditions of Lemma|42] Then, VM± € Z : M :- M_ + M+ > 2, V/Ts < M - I, the reconstructing polynomial on 
Si.M_.M+ (Proposition l2.2l l can be represented, almost everywhere, by combination of the reconstructing polynomials 
on the Ks-level substencils (Definition ! L2I ) of Si M_,M^, 



PRuM-mMi + ^Ax; X;, Ax; /) = ^ o-R^^M^^M^^K^kM) PRi,M--h,M+-K,+k(xi + ^Ax; Xi, Ax; f) 



V^eR\5R,,M_,M„;f, Vx,- e R VAx € R>o V/ : R ^ R 
where the rational weight-functions crR^_M_.M+.K,,kSS) '^''^ defined recursively by 

(^Ri,M-,M+-M-+k,M{^) 



O'R^,M_,M+,K^,kS0 



= 1 



Q'R,,M_-<:s,M+-l+its,-M-+fcM(^) 
mm(K^-l,k^) 

^ CRi,M-,M+,K,-i,eX^) '^Ri.M-~e„M+-(K,-r)+(„i,h--{X^) Ks>2 

{ 4=max(0A-l) 

V/Ts e {I,-- - ,M- 1) 



and satisfy the consistency condition 

its=0 

The set of poles of the rational weight-functions SR^^M-,Mt,Ki ( I39al l satisfies 

{24j 

<Sr,,M_,M^,1 : = {^ 6 R : Q'fi,,M_-l,M+.M+(^) = 0) = {^«,,M_-l,M+,M+.n; « 6 {"A^- + 1, ' ' ' , M+ - 1 



4=0 4=0 



Ls=0 4=0 



V/Ts e {I,-- - ,M - 1 



(39a) 



(39b) 



(39c) 



(39d) 



(39e) 



24 



Proof. The case - I follows from Lemma 14^21 with the set of isolated singular points >Sr,_m_.m+,i defined by 
( |39d| ). because of ( I32el ). Since the (K^ = l)-level subdivision rule is established, the conditions of |12, Lemma 2.1] 
are satisfied, proving (l39aH39cT i. and the recursive definition ( I39el i of the set of isolated singular points .Sr,,m_.m+,a:,- 
The c relation is used in (I39el i for > I, because there may be pole cancellation by the multiplications in (I39bb . □ 

Corollary 4.6 (Representation of the fundamental polynomials Q'R,,M-,M+,f(^) ( II lab ). Assume the conditions of Pm^o- 
sition 143] Then, the fundamental polynomials of Lagrange reconstruction ur^^m-M+A^) dl lab on Si^M-,M+ (Defi- 
nition II. lb , can be represented by a weighted combination of the basis (Proposition 13.3b Lagrange reconstructing 
polynomials on the Kg-level substencils (Definition ll.2l ) of Si^M_,M+ os 



fc=max(0/+A's-M+) 



(40) 



where the weight-functions a-R^ M-,M+,K^,kX^ '^'"^ defined by ( I39bb in Proposition 14.51 and the set of isolated singular 
points Sf>j M-,M+,K^ 

fcv (i39dll39el i. 



cip a; 
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Proof. Rewrite i 

uM-MMi + ^^x; Xi, Ax; f) ^ crRuM-M^,K„kX^) PR,,M--k„M^-K,+kAxi + ^Ajc; Xj, Ax; f) 

M+-K^+k, 



k,=0 

I—, 

IP 

L=0 



0-Ri,M_,M^,K„kS€) ^ [a;Rj^M_-k,,M^-K,+k,A^) f(Xi + £Ax)^ 



t=-M_+k. 



^ CRx,M-,M+.K^,kSS) 0:R\,M--k^,Mt-K,+k„e{^) 

Us=max(0/+A's-M+) 

a«i,M_,M+/(^) 



f(Xi+mx) (41) 



proving (|40ll by (l8a] i. 

Proposition 4.7 (Uniqueness of weight-functions o-R^^M-,M+,K„kM) ( I39bb ). Assume the conditions of Proposition 

The functions o-R,^M-,M+,K„kSO satisfying (I39ab are unique. 



□ 



Proof. We have proved by construction (|39bb the existence of weight-functions crR^ M_ M+.K„kSO satisfying ( |39a| 
I39cb . Uniqueness for the case /Ts = 1 was proved in Lemma W2\ Notice first that Corollary l4.6l does not require 
the validity of the particular expression (I39bb of the weight-functions, and is therefore valid for any set of weight- 
functions o-R,^M-,M+,K„kXO satisfying (I39ab . To prove therefore uniqueness we can use (l40b . which can be explicitly 



< k, <Ks M 
+ k,<e<M+-K, + k, ^ 



0<h<K, 

-M_ <e <M+ 

-M- +k, < C < - Ks + h 



0<ks<K, 
-M_ <f<M+ 
k,<t + M- 

e-M^ + K,<k, 



-M_ <e<M+ 

Q<k,<Ks 
e-M++Ks<ks<e + M- 
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written as 



V/Ts > 1 or^^m-.+m^-mA^) =o-Ri,M-,M^,K„oiO (^Ri,M-,m^-k„-mA^) (42a) 

V/Ts > 1 Q'B,,M_,+M+,-M_ + l(^)=0"fii,M_,M+,A:s,o(f) ^Ri ,M-,M+-A's,-M_ + 1 (f) 

+o"«,,M_,M+,is:s,i(f) Q'Ri,M_+i,M+-A:s+i,-M_+i(f) (42b) 

V/Ts > 2 a'fi,,M_,+M+,-M_+2(^)=0"«,,M_,M+,A:s,o(^) aRi,M_,M+-/s:s,-M_+2(f) 

+0"Ri,M_,M+,A's,l(^) a^«i,M- + l,M+-A's + l,-M_+2(^) 

+0"j?i,M_,M+,A's,2(^) 0'Ri,M-+2,M+-A:s+2,-M_+2(^) (42c) 



min(A:s/+M_) 
*:s=max(0,f+A:s-M+) 



V/Ts > 2 Q'R,,M_,+M+,+M+-l(^)=CrRi,M_,M+,A:s,A:s-2(^) Q'Ri,M_-A:s+2,M+-2,+M+-2(^) 

+o"Ri,m_,m+,a:s,a:s-i(^) Q'Ri,m_-a:s+i,m+-i,+m+-2(^) 

+o"Ri,a/_,m+,a:s,^^s(^) Q^Ri,M_-/f„M+,+M+-2(^) (42e) 

V/Ts > 1 ffR,,M_,+M+,+M+-l(^)=CrR,,M_,M+,A:s,A:s-l(^) 0;r^^m_-K^ + 1,M+-1,+M^-i(0 

+0'R,,M-,Mt,K,,KS^) 0!R],M--K^,M+,+Mt-l(0 (42f) 
V/Ts > 1 Q'R,,M-,+M+,+M+(^) =0"Ri,M_,M+,A's,/!:,(^) Q'Ri,M_-/!:„M+,+M+(^) (42g) 

Starting with (I42al |42g) we immediately prove uniqueness of o-R^^M^ M^.K^oiO cr^|,M_,M+,;i',,A:s(^), by contradiction 
because of ( |19b| l. Having proved uniqueness of crs,,M_,M+,;r„o(^), ( I42bl ) proves uniqueness of crsi,M_,M+,/r,,i(^), by 
contradiction because of (I19bb . In exactly the same way, having proved uniqueness of o-r^^m-,m+,k^,kSO' (I42fb proves 
uniqueness of o-R^ M-,Mt,K:„Ks-iiO- Continuing the procedure until reaching cr^j m_ m+ /fs [^](^) (f^'" increasing k^. 



starting from (I42al i) and cr^^ ^ ^ L'tJ*^^^ (for decreasing A;s, starting from (42gi), completes the proof of uniqueness. 

□ 

Corollary 4.8 (Weight-functions and approximation-errors). Assume the conditions of Proposition 14.51 Then, pro- 
vided that the reconstruction pair (Definition ! 1.3b of fix), h(x) :— [^(1;Aa)(/)](-^)> is sufficiently smooth Vjc e [x,_jv/_ 
- \ ■■ . x,+jv/+ + 5 Ax], for the expansions (I8bll8cl ) of the approximation error to hold, the weight-functions crR^^M-,M^,K^,k, (0 
(I39bl ) and the approximation-error polynomials fiR^M-,M^,n(0 ( II Ibb and AR^M-,M+,n(^) (lUcb . satisfy 



Yj CrRl,M-,M,,K,,kM) ^RuM^-k,,M,-KM.,iM) =0 V« E {M - /T^ + 1 , M} (43a) 

^ 0-RiM-,M+,K„kM) ^RuM--k„M+-K,+k„n(^) -'^RiM-M^AS) Sn> M + I (43b) 

^ 0-^,,M_,M+,/r,A(^) f^RuM^-k,.M,-K,+k,.n{0 =0 Vn E {M - + 1 , M} (43c) 

^ 0-R,,M_,M^,K„kM) mx,M_-k„M^-K,+k„n{^) =I^Ri,M_,M^,n(^) V« > M + 1 (43d) 

where the set of isolated singular points SR^^M_,M^,K^ is defined by ( |39dll39eT i. 
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Proof. The proof is quite obvious by replacing pr|_m_,m+(-^i+^Ax; x,-, Ax;/) and /?fi,,M_-i:s,M+-A:s+it,(-^i+^Ax; jc,-, Ajc; /) 
in ( |39a| ) by either (ISbT i or (|8c), yielding 

P«,,M_,M,(x,- + ^Ax; XI, Ax; f) ^h{x, + ^Ax) + ^ ^j?,,m_.m,.„(^) Ax" + ^Ax) + 0(Ax^^'^') 

n=M+l 



/z(x; + ^Ax) 



= 1 ( I39cl ) 

n=M-Ks + \ Us=0 

+6>(Ax^"+') 



Ax" h^"\xi + ^Ax) 



(44a) 



Pr„M-,mMi + ^Ax; Xi, Ax; /) ^h(xi + ^Ax) + ^ //«„m_.m,.«© Ax" /"^x,) + 0(Ax^"+') 



139^1 



^ a■R^,M-M^,K„kS^) PR^.M--k„M^-K,+K{xi + ^Ax; X,-, Ax; f) 



( JsTs 



Vt=0 



/!(X; + ^Ax) 



= 1 (I39cl) 

A'tj ( K, 

+ ^ 0"R,,M_,M+,A:„it,(^) I^Ri,M--h,Mt-Ks+ks,n(^) 

n=M-K, + l Us=0 

+0(Ax^"+i) 



Ax" f">(xd 



(44b) 



€ M \ 5jj„M_,M„jf Vx; E R VAx€M>o V/7eC^-+i(R) / := 
which prove ( l43T l by identification of coefficients of Ax". □ 
Example 4.9 (Rational weight-functions o'ii^^M_,Mi,Ki,k, 

(^) ( |39bl )). The stencil S;,3,3 (M_ = 3, M+ = 3, M M_ + 
M+ - 6) is symmetric around ^ = (Fig. [3]). The {Kf^ = j^yj = 3)-level subdivision (Definition II. 2l i is the highest 
level of subdivision for which all of the substencils contain either point / or point / + 1 (Fig. [3]). The rational weight- 
functions a-R. 3 3 3 (^) (fes £ {0, • ■ • , 3}) are all > in the interval / i around point ^ = +4 (Fig. |5]l. Because 

' ' ' ' ' Csj (^),3,3,3 ' ^ 

of the symmetry of the stencil 8, 3 3 around ^ = (Fig.|5]l, we also have cr^j 3 3 3 jt^(-i) > Vfcj e {0, • ■ ■ ,3}. The 
stencil S/3,4 (M_ = 3, M+ = 4, M := M_ + M+ = 7) is symmetric ai'ound ^ = 5 (Fig.©. The (K^ = ] = 4)- 
level subdivision (Definition 11.21 ) is the highest level of subdivision for which all of the substencils contain either 
point / or point / + l(Fig.|4]i. The rational weight-functions cr«| 3 44i^(f) {k^ e {0, ■ • ■ ,4)) are all > in the interval 

/ I around point ^ - +\ (Fig.|6ll. The stencil 8,3,4 not being symmetric around ^ = (Fig.|6]), positivity of the 

2)3,4,4 

weight-functions does not hold around ^ - where ctrj 3 444(-i) = -tTq, by direct computation using ( |39bt . The 
conditions of positivity of weight-functions at ^ = which is important in the development of WENO schemes 
are studied below ( il4.4l i. □ 
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b 




Ci?i,3,3,3,o(0 
C_Ri,3,3,3,l(0 
ffli,3,3,3,2(0 
O'i?i,3,3,3,3(0 
-^CR,(i),3,3,3 



Figure 5: Rational weiglit-functions o■|{^^M-,M+,Ks,k& I39bt for tlie (Kg = 3)-level subdivition (Definition II. 2t of tlie stencil S/3 3 (Fig.Q, and 



interval of convexity of the weight-functions around i - 



(Theorem|4l4). 



-a 



b 



1.5 - 



0.5 - 







-0.5 - 




fi?i,3,4,4,o(0 

o'fli ,3,4,4,1(0 

Ci?! ,3,4,4,2(0 
O'fli,3,4,4,3(0 
17^1,3,4,4,4(0 
4Rj(i),3,4,4 



Figure 6: Rational weight-functions o"/;, .m_,m+.a:s,*s(^) I39bl for the (Kg = 4)-level subdivition (Definition 1 1.2t of the stencil S,,3 4 (Fig.[2j, and 
interval of convexity of the weight-functions around (' - 



i / 1 (TheoremBj4). 



4.4. Convexity 

The nonlinear modification of the optimal (linear) weights in WENO schemes 10 [gl] is more straightforward when 
the combination dSal i is convex [120]. 

Remark 4.10 (Consistency, positivity and convexity). As can be seen by (l44all44bT i. condition ( I39cl ) ensures the con- 
sistency of the representation (I39ab as an approximation of h{x) =: [/?(1;Aa)(/)](-^) (Definition [OJ, and is therefore 
called the consistency condition of the representation ( |39a| ). Obviously, when at a fixed ^ e M all of the A's-level- 
subdivision weight- functions are > then, because of ( |39c| l, they must take values € [0, 1] (proof by contradiction) 

[(rR,.M-.MM(^) > V^s e {0, ■ ■ ■ , M [O < ctr.^m^^mmXO < 1 Vfe^ e {0, ■ ■ ■ , K,}] (45) 

Hence, positivity of the weight-functions at a fixed ^ e M ensures, by the consistency condition (I39cl l. that, locally, 
the representation (I39ab is convex. □ 
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In the early WENO papers 111 \m convexity of the combination ( fSal i had been postulated, and verified by direct 
determination of the coefficients at f = i nil. ShuH showed examples of combinations of choices of the stencil 
S,,M_,M+ (Definition ll.il ). of the level of subdivision ("Definition ll .2l i. and of the location ^ e M, for which convexity 
of ( fSal i is lost, and this appeared as a practical problem, not only in 2-D and 3-D unstructured grids ll2Cil . but also in the 
development of centered (central) WENO schemes (l^. For this reason the intervals of convexity were investigated 
numerically | 



a nil 



The analytical results obtained in the present work, in particular the recursive analytical expression of the weight- 
functions o-ji^ M-,M+,Ks.kAO (Proposition 14.5b and the factorization of the fundamental functions of Lagrange recon- 
struction Q'j?,,M_,M+,f(^) (Proposition 13.7b . can be used to study convexity intervals for arbitrary values of [M-t,Ks,], 
as was recently done for the Lagrange interpolating polynomial lfl2l. Proposition 3.2]. In |7i Result 6.1, p. 300] we 
had conjectured that for any choice of [M-t,Ks] for which all of the substencils Si^M_-k„Mt-K,+k, i^s 6 {0, ■ ■ ■ ,Ks}) 
contain either point / or point / + 1 (or both), convexity was observed at ^ = |. We provide here a formal proof of this 
conjecture, and give an estimate of the interval of convexity around ^ - \- 

Lemma 4.11 (Positive subdivision). Consider the subdivision level ^ 1 ofSi^M^M^ (Definition ll. 2b . Ijf 



- M_ < < 1 < M+ 
1<K,< min(M„ + 1,M+) 

then all substencils contain either point i or point i + 1 

HM^-hM,-K,+h n {/, / -H 1 } e[0,--- 

More precisely 



(46a) 
(46b) 



(46c) 



(I46all46bb 



[i, i+l] c S;,M_-t,M+-/f,+t Vfes e { 1 , ■ ■ • , - 1 ) 



(46d) 



A subdivision (Definition ll. 2b satisfying (I46db will be called a positive subdivision ^ Result 6.1, p. 300]. 



Proof. First notice that if all substencils contain either point / or point / + 1 (I46cb then so does the entire stencil 
Si,M-,M+ ^ Uf^o ^i,M--K,M+-KMt ■ Taking into account that by hypothesis > 1 , in the condition {/, / + 1 ) n S,;m_.m+ 9^ 
0, implies that S;_m_,m+ must contain both points / and / + 1 (prool^by contradiction taking into account > 1). The 
condition that both points {/, i + 1} must be contained in the big stencil S,,m_.m+ yields 



(2d) 



/,/ -H 1} n S,;m_,M+ + 

i,i + 1} <z [i - M ■ ■ ■ ,i + M+ 



(frjlD 



i-M-<i<i+\<i + M+ 



- M_ < < 1 < M+ 



(47a) 



proving that (146 ab is a necessary condition for the validity of (I46cb . Combining (I46cl 147 at implies (prool^by con- 
tradiction) that / must belong to the leftmost substencil (^s = 0) and ; + 1 must belong to the rightmost substencil 



^"Since by )46ct each of the substencils S,_M--/ts,Af+-A's+fe & - ■ ■ > i\) has a non-empty intersection with [i, i + 1), so does their union 
S,,M_,M+ {2dJ, ie \i, i + 1) n S,_m__m^ 0. Obviously the conditions {{i + M+ < ; < i + 1) V (;' < i + 1 < i - M_)) => {i, i + 1} n S,_m_,m+ = 
are a contradiction, implying that their negation is true, ie we must have ({i + M+ > i) A (i + 1 > i - It turns out that the inequalities in 

((i + M+ > ;•) A (i + 1 > I - M_)) must be strict. Assuming / + 1 — / — M_ i < i + \ < / — M_ + 1 {/, / + 1)0 ?ti^M--iM+-K^+^ ~ ^ 
contradicts l l46c> , implying / + 1 > i - M_ => M_ > -1 => M_ > 0. Assuming i = i + M+ =^ i+l > i > i + M+ - + (Ks - 1) =^ 

{i, i + I] n SiM-+{Ks-l),Mi--Ks+(Ks-l) = contradicts I46ct . because by hypothesis Kg > 1, implying i < i + M+ ==> M+ > => M+ > 1. 

' ' |473 , , BTa 

Assuming / i S,_m-,m+-a"s ' > i + M+ - => {!,('+ 1) n Sii M-,M+-Ki = <ZS contradicts (|46c) for kg = 0. Assuming i + I t SiM-+K|„^^-^ => 

i + \ < i - M- + Ks => {;',!+ 1 1 n S,j^_+Kj.m+ = contradicts )46cl for = >\. 
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(I46cll47al i 



/ + 1 e SiM_+K,.Mt 

>K, < min(M_ + 1,M+) 



i-M- < i < i + M+- 

i- + Ks < / + 1 < / + M+ 



(47b) 



proving that (I46bb is also a necessary condition for the vaHdity of ( I46cl ). To complete the proof it suffices to show that 
(146 al I46b] i are not only necessary but also sufficient conditions for (I46cb . We have 



( I46all46b] l 



M+ > 1 

-M_ < 

-M_ + < 1 

M+-K, > 

-M_ 

- M_ + < 



i- M- < i < i + 1 < i + M+ 

i- + < i+l 

i < i + M+ - Ks 



i-M^+K, < 

i - + 

i-M^+K, < 



< i <i + I < 

< i < 
i+l 

i < 

i+l < 

< i <i+l < 

< i < 

< i < i+l < 
i+l 



+ M+ 
+ M+ ■ 

+ M+ ■ 



+ M+ 
+ M+ 
+ M+ 
< / + M+ 



\lhe{Q,--- ,K,-l} 
K,+k, "ik.eil,--- ,K,] 



K, + k, "^k, e {I,--- ,K,- 1} 



completing the proof, the last conditions in (147 cl i being exactly ( I46dl 



(47c) 



□ 



Corollary 4.12 {{Ks = l)-level positively subdivisible stencils). Assume that M± e Z : M := M_ + M+ > 2 defining 
the stencil Si^M-.M+ (Definition ll.il) satisfy -M_ < < 1 < M+ ( I46al l. Then the {K^ - l)-level subdivision of the 
stencil S,_m_,m+ (Definition ll.2b is a positive subdivision (Lemma K.llb . 

Proof. By ( I46ab we have t hat ( (M_ > 0) A (M+ > 1)) ^ min(M_ + l,M+) > 1, so that - I < min(M_ + 1, M+). 
Hence the conditions ( 146 al I46bb are satisfied, so that, by Lemma |4. Ill the {K^ = l)-level subdivision of a stencil 
satisfying ( 146 ab is a a positive subdivision. □ 



Lemma 4.13 (Convexity in the neighborhood of / + j for (Ks = l)-level subdivision). Assume that M± e Z : M :- 
M_ + M+ > 2 ( l2a] l defining the stencil S,;m_,m+ (Definition ll.il ) satisfy -M_ < < 1 < M+ ( I46ab . Then the rational 
weight-functions crs,,M_,M+,i,o(^) ( I32ab flw^/ crRi,M_.A/+,i,i(^) ( 1 3 2bl ) for f/ie representation of the Lagrange reconstructing 
polynomial by the Lagrange reconstructing polynomials of the (Ks — l)-level substencils of Sj M^.Mt (Lemma 14.2b 
satisfy 



< 0"R,,M_,M+,l,fe(^) < 1 



C«,(i),M_,M_ 



+,1 " (^C„,(i),M_,M+,r^C„,(|),M_,M+,l) ■ 



(48a) 



Vfcs e {0, 1} 



where the limits of the convexity interval around ^ ■ 



5> ^Cfl (-_) M- M+ 1 ^ 5 of length > 0, are defined by 



'Cfl,(i),M_,M+,l 



^C„,(i),M_,M+,l 



max(^R,_M_,M+,+M+,0,^^i,M_,M+,-M_,Oi^«i,M-,M^_-l,-M_,o) -W- >0 



M+ = 1 
M , > 1 



(48b) 
(48c) 



where ^Rt.M-.M+,e,n (n e {-M_, • ■ • ,M+} \ are the M real roots (Proposition [33]l of the fundamental polynomial of 
Lagrange reconstruction aRt^M-,M+AO ( 1 1 1 
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Proof. By hypothesis, the stencil S,;m_,m+ satisfies the conditions of Corollary 14. 121 implying that the (K^ - l)-level 
subdivision of S,;m_,m+ is a positive subdivision, satisfying the conditions of Lemma l4.11l and we have by ( I46dl l 



{/,/+l)C S,_M_,M+ 



(49a) 
(49b) 
(49c) 



By Lemma WIT] the rational weight-functions cr«, m_,m+,i,o(^) i32a\ and 0-^, m_,m+,i,i(^) (I32bl l can be expressed in 
terms of the fundamental polynomials of Lagrange reconstruction (Proposition l2.2l l Q'r,,m_,m+,-m_(^), Q'«,,m_,m+,+m+(^), 

Q'«i,M-,M+-i,-M-(^), andaS[_M_-i,M+,+M+(^)- Notice that, because of the identities of Proposition l3 . 1 Ol we have aR,^M--i,M+, 
Q'fi,_M_,M+-i,-M_(^)- By Proposition |33] all of the roots of the fundamental polynomials of Lagrange recon- 



struction are real, and therefore the factorization of Proposition 13.71 applies. Applying the factorization ( I25ab . and 
taking into account (I49al I49bl ). which were shown in Lemma 14111 to be direct consequences of (I46ab . we have 



M 



Ml 



M! 



<0 Vf <ft, 



n=l 

M+ 



,M_,M+,-M_,n) 
,M-,M+,-M-,;i) 



M_ = 



M_ > 



(49d) 



>0 Vf>fR| jn_ M+.-M-.O 



<0 Vf<f„, 



Q:RuM_,M+,+M^(g) - — I I - g«,,M_,M+,+M+,n; 



l4fal(-l) 



2M+ 



M! 



>0 V^>^^. M_,Af+.+M+.0 



n=-M- 




M+-1 



W - ^Ri,M_,M+,+M+,n) ]~~[ (^ - ^R,,M_,M+,+M+,n) 



M+ = 1 
M+ > 1 



(49e) 



>0 '^^>^Ri,M-~M+,+M+.0 



n=l 



<0 V^<^flj^^_^^^_|_^+^^_|_j 



M-1 ^+- ' 



(M-iy 



«=-M_ + l 



l4fal(-l) 



M-1 



(M-1)! 



<0 V^<^Rj 



^+-i,-M_,«) M_ = 



>0 'i^>^R-^,M-.M+-l-M-.Q 



,M_,M+-1,-M_ 







«=1 



]~~[ (^ - ^Ri,M_,M+-l,-M_,«) ]~~[ (^ - ^«,,M_,M+-1,-M_,n) 



>0 V^>^fij ,M_,M+-1,-M_,0 



<0 V^<^/;j 



M+ = 1 

M_ M+ - 1 
(49f) 



where in ( I49fb we only need to distinguish 3 cases because the constraint M :- M_ + M+ > 2 ( l2a] l implies 
that we cannot have simultaneously M_ - and M+ = I. Since by Proposition 13.101 aR,,M--\,M+M+(^) 

Q-R,_M_,M,-i -M_(f) ^ ^RuM-,M+-i-M-,ii = ^Ri,M_-i,M„M„« Vn € {-M_ -I- 1 , ■ ■ • , M+ - 1 }, defining the limits 
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of the open convexity interval around ^ - \, Iqr (■!■) M- m+ i < l48al i. by ( I48bll48c] i. we have 

sign(a«„M_,M„-M_(^))'^(-l)^"^^'i'(-l)''- ^ IcR,iM,M_,M.,i (49g) 

sign(a«„M_.M,-i,-M_©)'^(-l)^-'"^-' 'i'(-l)''- 6 4^_(,)_^ ,^^,i (49i) 



(_iy™ . e /c„,(i),M_,M„i (49j) 



sign(aR,,M_-i,M+,+jM+(f)) ="(-!)*' 'sign(Q'R,,M_,M+-i,-M_(f)) 

l49fl 



whence 



I32ill49g||49il 

cr«„M_,M.,i,o© ^ V^e/^^_^,^^^^^ (49k) 

I32bll49hl l49il 

cr«„M_,M.u(f) >^0 Vfe/^^_^,^^^^^ (491) 

Because of the consistency condition ( I32cl i. positivity of the weight-functions implies convexity (Remark [4.10l i. so that 
(|49k[ |491t prove ( 148 ab . Notice that, by Proposition |33] (I24ab implies that ^fi,,M_,M+/#o,o < 5 and ^R,,M_,Mt,f#i,i > 5, 
VM± € Z : M := M_ + M+ > 2, satisfying the conditions of Corollarv 14 121 so that the length of I^^ (i)M- m+ 1 
( I48ay48c] i is > 0. ' ' ' ' □ 

Theorem 4.14 (Convexity of a positive subdivision in the neighborhood of / + |). Assume that the subdivision level 
Ks > I ofSiM_,M+ (Definition ll.2l ) satisfies the conditions of Lemma l4. 1 \\ ( positive subdivision), viz 

M M_ + M+ > 2 ^ 
- M_ < < 1 < M+ ( l46al i 
1<K,< min(M_ + 1, M+) ( l46b] i 

implying (Lemma l4.11l ) f/zaf a// substencils contain either point i or point i + \. Define the interval 

^Cgj(^),M-,M+,A:s {~\ ^CKj(i),M_-4,M+-Ls+4,l (50a) 
L,=0 4=0 

recursively using convexity intervals /^^ (i)M_ m+ 1 (l48alj48cT l of {K^ — l)-level positive subdivisions (Lemma 14. 13b . 
Then the rational weight-functions a-R^ M-,Mt,K^,kS^) ( I39bl) satisfy 

V/Ts 6 {I,--- ,M- 1} (50b) 



< 0-R,,M_,M^,K„kM^ < 1 



implying convexity of the combination (|39ab . 

Proof. The vaUdity of (I50all50bb for A's = 1 was proven in Lemma I4TT3] Assume min(M_ + 1, M+) > 2 so that the 
(A's = 2)-level subdivision be a positive subdivision (Lemma [4.11b . Then, by (I39bb . we have 

O'R,,M_,M+,2,0(^) ' = 'o'«i,M_,M+,l,()(^) 0"R|,M_,M+-l,l,o(^) (51a) 

+0"«,,M_,M+,l,l(^) 0"«,,M_-l,M+,l,o(f) (51b) 
0"fii,M_,M+,2,2(f) '^Crs,_M_,M+,l,l(^) 0"«, ,M_-1,M+,1,1 (f) (51c) 
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Having assumed that the (Ks - 2)-level subdivision is a positive subdivision (Lemma [4.11l ). we have 



min(M_ + 1,M+) > 2 



M_ + 1 > 2 :S -(M_ - 1) < < 1 < M+ 
M+>2M -M_ < < 1 < (M+ - 1) 



(5 Id) 



implying that the 1-level subdivisions of the stencils S,,m_-i,m+ and S,;m_,m+-i are positive (Corollary 14. 12) . Therefore 
all of the 1-level weight-functions on the RHS of (I51al451c] i are positive in the neighborhood of ^ = j, because of 
Lemma l4.13l and we have 

< O-R,,M_,M„2,0© € (/ 

< CTRj^ ^^ 2,2 



1 n/ 1 

Ck, (5),M_,M+,1 C«, ( 5 ),M_,M+- 1,1/ 

1 n/ 1 n/ 1 1 

C«,(i),M_,M+,l C„,(i),M_,Mt-l,l C„,(^),M_-1,M^,1/' 



n/ 



Ck,(2),M_,M+,1 C„,(2),M_-I,M. 



Defining 



^Cfi,(l),M_,M+,2 I 11 I ^( 



1 



C«, (4),M--4,M+-Ls+4,1 ^c„,(i),M_,M+,l ^Cs,(i),M_,M^_-l,l ^Cs, (i),M_-l,M+,l 



(51e) 
(51f) 
(51g) 

(51h) 



we have that all of the 3 (K^ = 2)-level weight-functions (I51el -51gi are simultaneously positive e /^^ (i)M_ m+ 2' 
which (Remark k.lOb . because of the consistency condition (I39bl ). proves ( I50al ISObT i for - 2. 

It is straightforward to complete the proof by induction. Since we have already proved (I50al 150b) for A's = 2, 
assume - I >2 <=> > 3. By Proposition |4~: 



minCATs-l.jts) 



)39bl V-i 

4=max(0,jt5-l) 

Assume that (I50all50bl i are valid for - 1 > 2 

^^^■^C«j(i),M_,M„/f,-l " n n^C«,(i),M_-4,M,-L.+4,l 



(51i) 



< 0-R^^M-,M+,K,-l,eSS) < 1 



Z,=0 4=0 



(51j) 



t V4 e {0,/:^ - 1) 



and that A's satisfies (146b I ), ie 

min(M_ + 1,M+) > K, 



M_ + 1 > /Ts 



IS} 



-(M_ - ^s) < < 1 < (M+ - (/Ts - 1) + 4) 

V^-s e {O,--- ,/:s - 1) (51k) 

so that the substencils S,;m_-4,m+-(a:s-1)+4 satisfy the conditions of Corollarv l4.12l implying by Lemma l4.13l that 

n \Ji- i ^"^s^lO, 1} ,,,,, 

< cr^„M_-CM,-(^:,-i)+4.i.,«,(^) e ^^,^(1,,^ | y^^ e {0, . . . _ 1} (511) 

Combining (ISTTl [51j] ISTTli yields 



< CrR,.M-.M^.K,.kX^) 



'cS|(l),M_-4,M+-(X,-l)+4,l 



h'c 



V4=0 



i V^se{0,/:s} 

which (Remark l4.10b . because of the consistency condition (I39bl ). proves dSOal I50bb V/^s satisfying ( I46bl) 
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(51m) 



□ 



1.5 



1 - 



0.5 



- 



-0.5 



1 



M = 2k ; keN 
M = 2k + l; keN 



-| 1 1 1 1 — TT 



-| 1 1 1 1 — TT 



C«,(i),[f J,M-[f J,[f] 





_l I I I I 



-I I I I I l__L 



10 



100 



Figure 7: Interval of convexity /^^ (j) [^J ^-[7] f "1 ' 5 (Theorem 14. 141 . of the maximum positive subdivision-level (Lemma l4.1 U 

Ks = IyI' °^ usual WENO stencils S^. 1 m 1 1 m i (Definition ll.U . as a function of stencil width M (in logscale). 



Example 4.15 (Convexity around / + ^ of usual WENO discretizations). The usual WENO discretizations for the nu- 
merical approximation of f'(x) lllllil use 10, p. 298] the (K, = ])-level subdivision (Definition O of the 
general family of stencils S, |^mj ^_|^mj, which (Definition ll.il ) contains M + 1 points. If M = 2k (k e N>o) is even, 
then the stencil S, j m_|^m j - S,,^,^ is symmetric around point /, and upwind-biased with respect to the cell-interface 
/ + I (eg S,-,3,3; Fig.|3]l, corresponding to the family of WENO(2r - 1) (r := A; -1- 1) upwind-biased schemes I9I [loll . 
If M = 2A; -I- 1 (A; € N>o) is odd, then the stencil S,- j m-[-^J - Si^k,k+\ is symmetric around the cell-interface ; -1- 5, 
and downwind-biased with respect to the point / (eg S, 3 4; Fig.|6]l, corresponding to centered (central) WENO schemes 
ifTsil . For the family of stencils S,- j ^_|^mj, we have 



M. =M- 



min(M_ -1- 1 , M+) = 



M 



min(A; + l,k ) - k 


VM 


^2k 


AceNx) 1 


M- 


mm(k +l,k+l)-k+l 


VM 


= 2k+\ 


ken^Q ] 


2" 



VM E N> 



(52a) 
(52b) 
(52c) 



so that, by Lemma l4. Ill TT^ = min([^-y|-i4,M-|^Yj) ^ |^yj corresponds to the maximum level of positive subdivision 
of S, J M-[^ J ■ Therefore, Theorem l4. 14l applies. and there exists an interval of convexity /^^ (i) [- J m-[-!^J [^] < l50al l 
around^ = | (Fig.|7]l. 

Notice that the interval of convexity 7^,^ |^mj m-[^J \^■'\ (Fig-ID, is slightly larger for M - 2k (k e N>o) even, 
compared io M - 2k + \ (k e N>o) odd, and its length slightly decreases (quasi-logarithmically VM € {2, ■ ■ ■ ,22)) 
with increasing number of cells in the stencil, M (Fig. |7]i. For stencils with M - 2k (k e N>o) even, like 8,33 
(Fig.EJ, because of symmetry with respect to point it is s traig htforward to show that there is a symmetric interval 
of convexity around / - 5 (Fig.|5]l, as was also observed in fll. Tab. 3.2, p. 516]. On the contrary, for stencils with 
M - 2k+l (k e N>o) odd, like 8,3,4 (Fig.lH, it turns out that positivity of the weight-functions does not hold at ^ = -5 

□ 



(Example 14 . 9 1 Fi g . I&t . as was also observed in [ 1 1 , Tab. 3.5, p. 518]. 
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5. Conclusions 

In the present work, we studied analytically the representation of the Lagrange reconstructing polynomial by 
combination of substencils, and in particular the conditions under which this representation is convex, ie the weight- 
functions E [0, 1]. 

We first formalized several results on the fundamental polynomials of Lagrange reconstruction (Proposition 12.2) . 
o:r,,m-.m+AO < I1 lab . Each of the polynomials aR^^M-,M^,e{S) the reconstruction pair (Definition 11.3b of the corre- 
sponding fundamental function of Lagrange interpolation ai M_,M^,[{S) (Proposition |TlJ, and for this reason all of its 
M roots are real (Proposition |33), distant < | from the corresponding root of the fundamental function of Lagrange 
interpolation a/,M_,M+/(f) ( I12al ). This leads to a simple factorization of the fundamental polynomials of Lagrange 
reconstruction (Proposition [Jj)- 

The leading (9(Ax*^) term of the approximation error of the Lagrange reconstructing polynomials on 2 overlapping 
stencils shifted by 1 cell, {/ - M_, ■■■ ,i + M+ - 1} and {/ - M_ + \,--- ,i + M+}, is different (Proposition 13. 10b . and 
several identities hold between some of the fundamental polynomials on the 2 stencils. Based on these identities 
(Proposition ^ . 1 Ob . we show that there exist unique rational weight-functions combining the Lagrange reconstructing 
polynomials on {/ - M_, • ■ ■ ,i + M+ - 1 ) and {/ - M_ + I, - ■ ■ ,i + M+ ] into the Lagrange reconstructing polynomials 
on {i - M_, • • ■ ,i + M+] (Lemma [4.2b . this representation failing at the poles of the weight-functions, all of which are 
real and can be identified with roots of fundamental polynomials of Lagrange reconstruction. Having established this 
1 -level subdivision rule, the general recurrence relation for the weight-functions proven in Lemma 2.1] applies, 
and provides the analytical expression of the weight-functions for a general level of subdivision (Proposition l4.5b . and 
of the set of their poles, all of which are real. These weight-functions are unique (Proposition [47). 

Finally, we prove (Theorem l4.14b that for any K^-level subdivision of {/ - M_, ■ • ■ ,i + M+] into A's H- 1 substencils 
{i - M_ + kfi,- ■ ■ , M+ - + k^] {ks e |0, • ■ • , Ks}), iff each of the substencils contains either point / or point / + 1 
(positive subdivision; Lemma |4. lib , then there exists a neighborhood of ^ = ^ {x - Xi + |Ax), 7,^^ m+ k. ^ 5' 
whose limits can be explicitly defined by roots of fundamental polynomials of Lagrange reconstruction, where all 
of the weight-functions ctr^^m^m+.k^A^^^ > e 7;,^ {1)m_m^k^ implying because of the consistency relation 
2^*10 crs,_M_,A/+,A'sA(^) = 1 e K ( |39cb . that the representation of the Lagrange reconstructing polynomial by 
combination of substencils is convex e 7^^ m- m^k^ ^ 5- Theorem l4. 14| provides a formal proof of (and general 
conditions for) convexity in the neighborhood of^- |, which had always been conjectured, on the basis of numerical 
evidence, all along the development of WENO schemes lflll2il5il9i[3i [Tll[loll . 
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